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SECTION  I 


•  INTRODUCTION  AND  HISTORICAL  REVIEW 

1.1  INTRODUCTION 

There  has  beer,  an  increased  emphasis  directed  toward  refinements  in 
the  design  of  concrete  members.  For  refinements  to  be  possible,  a  know¬ 
ledge  of  the  behavior  of  concrete  under  combined  stresses  is  a  necessity. 
Past  research  on  concrete  behavior  under  combined  stresses  has  generally 
been  limited  to  the  study  of  the  relationships  between  the  applied 
stresses  at  failure.  Information  on  the  relationships  between  the  ap¬ 
plied  stresses  and  the  resulting  strains  was  limited.  None  of  the 
investigators  of  concrete  behavior  reported  any  attempts  to  mathemati¬ 
cally  relate  stresses  to  strains  for  the  combined  stress  conditions. 

This  contract  was  awarded  to  New  Mexico  State  University  to  perform 
an  experimental  investigation  on  the  behavior  of  plain  concrete  under 
various  biaxial  and  triaxial  loadings.  Some  of  the  terms  of  the  contract 
were  as  follows: 

(a)  Develop  the  necessary  techniques  and  perform  the  necessary 
tasks  to  record  the  strains  in  the  three  principal  directions 
of  loading. 

(b)  Construct  or  obtain  a  loading  device  which  is  capable  of  pro¬ 
ducing  a  maximum  compressive  stress  of  30,000  psi  and  a  tensile 
stress  to  failure. 

(c)  The  loading  de rice  will  be  capable  of  three-dimensional  loading. 

(d)  Concrete  with  a  nominal  unconfined  compressive  strength  of 
4000  psi  will  be  used  in  all  specimens. 
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(e)  The  applied  loads  will  range  from  a  tensile  load  causing  failure 
to  a  maximum  compressive  load  producing  a  stress  of  30,000  psi 
in  the  test  specimen. 

(f)  Develop  relations  which  will  enable  prediction  of  the  state  of 
stress  for  a  concrete  element  from  a  state  of  strain. 

(g)  Determine  a  failure  theory  that  is  consistent  with  observed 
failure  modes. 

(h)  Develop  a  computer  program  which  incorporates  the  failure  theory 
and  yields  a  set  of  stresses  for  a  given  set  of  strains. 

1.2  HISTORICAL  REVIEW 

Many  investigators  have  reported  experimental  and  theoretical 
studies  on  the  failure  characteristics  of  concrete  subjected  to  combined 
compressive  and  tensile  stresses.  Most  of  the  investigators  concentrated 
on  compressive  loadings  only.  All  of  the  publications  listed  in  the  List 
of  References,  with  the  exception  of  (19),  (20),  and  (24)*,  contain  in¬ 
formation  on  concrete  behavior  studies.  Generally  an  attempt  was  made  to 
associate  the  failure  stresses  to  an  existing  failure  theory  or  to  a 
modification  of  an  existing  failure  theory.  Association  of  the  failure 
stresses  to  an  existing  failure  theory  was  never  completely  successful. 

Failure  theories  that  have  been  considered  are  listed  below  and  are 
discussed  on  the  basis  of  their  agreement  with  observed  results  from 
tests  on  concrete  cubes. 

(a)  Mohr's  Failure  Theory  -  The  failure  mode  for  this  theory  is 


* 

Number  in  parenthesis  refers  to  List  of  References. 
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slippage  along  a  plane  inclined  with  respect  to  the  principal  stresses. 
The  observed  failure  mode  was  splitting  of  the  cubes  such  that  the  frac¬ 
tured  surfaces  were  essentially  perpendicular  to  the  direction  with  the 
lowest  applied  stress.  In  addition,  the  shear  on  a  slip  plane  is  assumed 
to  be  a  function  of  the  maximum  and  minimum  stresses  only  and  does  not 
consider  the  effect  of  the  intermediate  stress.  Test  results  indicate 
that  the  intermediate  stress  does  effect  the  failure  envelope. 

(b)  Octahedral  Shear  Theory  -  This  failure  theory  assumes  that 
failure  occurs  whenever  the  shear  stress  on  planes  whose  normals  possess 
equal  directional  cosines  with  respect  to  the  principal  stress  axes 
reaches  a  specific  value.  It  is  impossible  to  relate  the  octahedral 
shear  stress  to  a  failure  mode  and  it  has  not  adequately  described  the 
experimental  results, 

(c)  Distortion  Energy  Theory  -  This  theory  assumes  that  failure 
will  occur  whenever  the  energy  of  distortion  exceeds  a  constant  value. 

The  energy  of  distortion  is  the  difference  between  the  total  energy  in  an 
element  and  the  energy  due  to  a  hydrostatic  loading  (equal  stresses  in 
all  three  directions).  Computation  of  the  energy  of  distortion  is  not 
convenient  and  it  does  not  agree  with  experimental  results. 

(d)  Maximum  Tensile  Strain  Theory  -  Failure  occurs  whenever  the 
largest  tensile  strain  exceeds  a  constant  value.  The  failure  mode  of 
this  theory  conforms  to  that  observed  from  tests;  however,  there  is  not 
agreement  between  the  failure  stresses  predicted  by  this  theory  and 
experimental  results. 

(e)  Griffith's  Theory  -  This  theory  predicts  a  failure  mode  similar 
to  the  experimental  results.  It  considers  the  existence  of  microscopic 
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cracks  and  the  propagation  of  these  cracks  due  to  high  tensile  stresses 
generated  at  their  tips.  It  requires  an  analytical  description  of  the 
stresses  around  the  cracks;  hence,  its  use  is  limited  and  an  experimental 
verification  is  difficult. 

(f)  Maximum  Compressive  Distortion  Stress  Theory  -  This  theory 
assumes  that  failure  occurs  whenever  the  compressive  stress  which  produces 
distortion  of  an  element  exceeds  a  constant  value.  The  distortion 
stresses  are  the  difference  between  the  stresses  acting  on  an  element 

and  the  average  of  all  of  the  stresses  acting  on  the  element.  This 
theory  does  not  agree  with  experimental  results. 

(g)  Maximum  Tensile  Distortion  Stress  Theory  -  This  theory  assumes 
that  failure  occurs  whenever  the  tensile  stress  of  distortion  reaches  a 
limiting  value.  The  distortion  tensile  stress  is  obtained  in  the  same 
manner  as  mentioned  for  the  compressive  distortion  theory.  The  theoret¬ 
ical  failure  mode  agrees  with  the  experimental  failure  mode;  however,  the 
failure  stresses  predicted  by  this  theory  do  not  agree  with  the  experimen¬ 
tal  results. 

(h)  Modification  of  the  Maximum  Tensile  Distortion  Stress  Theory  - 
Zimmerman  (29)  proposed  an  empirical  modification  of  the  maximum  tensile 
distortion  stress  theory.  Combining  the  empirical  modification  and  the 
tensile  distortion  stress  theory  resulted  in  adequate  agreement  between 
the  failure  stresses  whenever  compared  to  experimental  results.  This 
agreement  was  limited  to  a  rather  limited  range  of  principal  stresses. 

A  review  of  the  existing  failure  theories  as  related  to  concrete 
behavior  leads  to  the  conclusion  that  they  cannot  be  used  directly  to 
describe  the  failure  characteristics  of  concrete. 
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A  review  of  publications  on  concrete  behavior  under  combined  stresses 
revealed  that  only  a  few  of  the  investigators  reported  the  results  of 
strain  measurements.  In  no  case  was  there  reported  any  attempt  to  obtain 
relationships  between  stresses  and  strains. 

1.3  SCOPE  OF  PRESENT  STUDY 

The  present  study  was  designed  to  determine  experimental  stresses 
and  corresponding  strains  for  combined  stress  conditions.  The  combined 
stress  conditions  included  both  compressive  and  tensile  stresses  as  well 
as  all  combinations  of  compressive  and  tensile  stresses  acting  together. 
This  also  includes  uniaxial,  biaxial,  and  triaxial  states  of  stress. 
Presently,  dynamic  codes  calculate  a  stress  matrix  from  a  given  strain 
matrix.  Therefore,  this  effort  will  consist  of  the  formulation  of  con¬ 
stitutive  relationships  such  that  a  set  of  stresses  can  be  determined  from 
a  given  set  of  strains.  In  experimental  considerations  the  strains  would 
be  measured  and  the  stresses  computed  from  the  measured  strains. 

1.4  STATE-OF-THE-ART  IN  CONCRETE  TESTING  AND  MODELING 

A  review  of  the  published  results  indicate  that  several  different 
approaches  have  been  taken  in  the  testing  of  concrete  subjected  to  multi- 
axial  loadings.  Multiaxicl  tests  were  conducted  using  hollow  and  solid 
cylinders,  cubes  and  square  slab-type  specimens. 

Hollow  cylinders  have  been  used  in  specimens  for  biaxial  compression, 
biaxial  tension-compression  and  biaxial  tension  tests.  For  biaxial  com¬ 
pression  tests  using  the  hollow  cylinders,  a  triaxial  state  of  stress  is 
suspected  and  the  reported  strengths  are  considered  to  be  too  high. 
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Solid  cylinders  have  been  used  as  specimens  for  biaxial  and  triaxial 
compression  tests.  This  type  of  specimen  has  a  disadvantage  since  it  is 
not  possible  to  independently  control  the  applied  stress  in  three  direc¬ 
tions  . 

Cubes  have  been  used  primarily  for  biaxial  and  triaxial  compression 
test  specimens.  When  this  type  specimen  was  used,  with  no  provision  made 
for  the  reduction  of  friction  between  the  specimen  and  the  loading  platen, 
the  reported  results  are  considered  to  be  too  high. 

Several  investigations  have  been  conducted  using  square  slab-type 
specimens.  Biaxial  compression,  biaxial  tension-compression  and  biaxial 
tension  type  tests  have  been  conducted  using  this  type  of  specimen.  This 
type  of  specimen  seems  to  develop  a  reasonably  unconfined  biaxial  stress 
state.  This  type  of  specimen,  however,  does  not  seem  to  be  suitable  for 
a  triaxial  test  specimen. 

There  are  very  few  published  results  in  connection  with  the  modeling 
of  concrete  behavior  for  concrete  subjected  to  multiaxial  loading. 

Anson  (23)  used  a  model  to  simulate  the  behavior  of  concrete.  His 
model  consisted  of  nodal  points  connected  by  bars  of  different  stiffnesses. 
The  nodal  points  are  located  at  the  points  of  an  octahedron.  Bars  of  a 
constant  stiffness  ratio  were  used.  The  model  was  used  to  simulate 
stress-strain  curves  of  concrete  in  uniaxial  tension  and  compression  and 
the  failure  envelope  of  concrete  subjected  to  biaxial  loading  and  triaxial 
loading.  The  author  concluded  that  his  model  should  only  be  used  for 
qualitative  examination  of  concrete  behavior. 

The  results  of  an  experimental  and  analytical  investigation,  of 
simplified  models  of  concrete  to  study  the  strength  and  deformational 
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behavior  of  plain  concrete  subjected  to  short  term  compressive  loads  was 
reported  by  Buyukoztush,  Nilson  and  Slate  (32).  The  study  was  limited  to 
uniaxial  and  biaxial  states  of  stress.  Thin  square  plates  of  mortar  with 
round  stone  inclusions  of  uniform  size  were  loaded  in-plane.  The  experi¬ 
mental  results  were  compared  to  a  finite  element  analysis.  The  model  used 
had  a  disadvantage  in  that  the  solution  involved  428  triangular  finite 
elements  to  model  the  behavior  of  a  5  x  5  inch  specimen. 

The  model  did  not  predict  the  nonlinear  load-deformation  curves  very 
accurately.  The  authors  acknowledged  that  the  two-dimensional  representa¬ 
tion  of  the  material  adopted  for  their  study  was  not  entirely  satisfactory 
even  for  the  two-dimensional  stress  state  considered. 
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SECTION  II 


MATERIALS  AND  MIXING  PROCEDURES 

2.1  CONCRETE  MATERIALS 

The  water  used  in  mixing  the  concrete  was  taken  from  the  water  system 
supplying  New  Mexico  State  University.  Type  I  cement  was  used  and  was  pur¬ 
chased  locally.  No  control  was  exercised  over  the  quality  of  the  water  or 
the  cement. 

The  fine  aggregate  (sand)  and  the  coarse  aggregate  (gravel)  were 
both  purchased  locally  and  from  the  same  supplier.  The  sand  was  taken 
from  the  supplier's  stock  and  delivered  to  the  laboratory  and  placed  in  a 
weather-protected  storage  bin. 

The  gravel  required  for  this  project  (maximum  size  of  aggregate  was 
limited  to  1/2  inch)  was  not  stocked  by  any  local  supplier;  therefore,  it 
had  to  be  specially  produced.  Approximately  five  cubic  yards  of  the  gravel 
were  processed  and  delivered  to  the  laboratory.  The  gravel  was  stored  out¬ 
side  the  laboratory  and  was  exposed  to  the  weather. 

The  results  of  sieve  analyses  on  the  sand  and  gravel  are  shown  in 
Figure  1.  The  variation  in  the  sieve  analyses  is  indicated  along  with 
the  average  of  all  the  sieve  analyses. 

All  results  from  sieve  analyses  of  the  sand  fell  within  the  ASTM 
grading  limits.  There  was  variation  from  test  to  test;  however,  this  is 
mostly  due  to  the  difficulty  in  obtaining  representative  samples  to  be 
used  in  the  sieve  analyses.  The  sand  was  removed  from  the  bottom  of  the 
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storage  bin.  Some  segregation  of  the  sand  sizes  no  doubt  occurred  at  the 

time  the  sand  was  dumped  into  the  storage  bin.  The  larger  particle  sizes 

roll  to  the  outside  of  the  sand  column  and  the  fines  tend  to  concentrate 

in  the  center  of  the  sand  column.  The  sand  sample  was  obtained  by  taking 

♦ 

two  or  three  shovelsful  of  the  sand  and  reducing  the  amount  with  the  aid 
of  a  sand  splitter.  The  sample  weights  ranged  from  935  to  1710  grams. 

The  Fineness  Modulus  of  the  sand  was  also  computed  for  each  sieve 
analysis.  The  average  Fineness  Modulus  was  2.69  and  the  range  was  2.55 
to  2.92.  The  entire  range  was  within  the  recommended  Fineness  Modulus 
limits  of  2.3  to  3.1.  The  Fineness  Moduli  for  the  various  samples  are 
listed  below. 


tmple 

Fineness  Modulus 

1 

2.65 

2 

2.55 

3 

2.77 

4 

2.66 

5 

2.74 

6 

2.59 

7 

2.74 

8 

2.63 

9 

2.68 

10 

2.92 

11 

2.64 

In  two  of  13  tests,  the  sieve  analysis  results  for  the  gravel  exceeded 
the  ASTM  grading  limits  for  the  No.  4  size.  The  average  of  all  the  test 
results  did  fall  within  the  ASTM  grading  limits.  It  was  intended  that  no 
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batch  falling  outside  the  ASTM  grading  limits  be  used.  Obtaining  repre¬ 
sentative  samples  of  the  gravel  for  sieve  analysis  was  far  more  difficult 
than  in  the  case  of  the  sand. 

Samples  of  gravel  to  be  used  in  the  sieve  analyses  were  obtained  by 
taking  two  or  three  shovelsful  from  a  batch  and  reducing  the  sample  size 
by  using  the  sand  splitter.  The  gravel  would  pass  through  the  sand  split¬ 
ter.  Gravel  sample  weights  ranged  from  840  to  1805  grams. 

The  gravel  contained  more  than  the  specified  maximum  percentage  of 
flat  and  elongated  particles.  Twenty  percent  maximum  was  originally 
specified,  however,  the  gravel  contained  approximately  twenty-five  per¬ 
cent  flat  and  elongated  particles.  The  gravel  consisted  of  crushed 
stone.  It  appears  that  the  flat  and  elongated  particles  resulted  from 
the  crushing  operation. 

2.2  CONCRETE  MIX  DESIGN  CRITERIA 

The  concrete  used  in  forming  the  cube  test  specimens  was  designed  on 
the  basis  of  a  28-day  strength  of  4000  psi  and  a  slump  of  3-4  inches.  The 
rather  large  slump  was  required  to  form  the  relatively  small  cube  test 
specimens.  Good  workability  was  also  Important.  One  side  of  the  cube 
specimens  had  to  be  manually  smoothed  with  a  trowel. 

2.3  TRIAL  CONCRETE  MIXES 

The  first  trial  concrete  mix  was  proportioned  on  the  basis  of  the 
absolute  volume  method  recommended  by  the  Portland  Cement  Association 
iPCA) .  The  material  quantities  for  the  first  trial  mix  was  estimated 
using  charts  and  graphs  published  by  the  PCA.  The  design  mix  was  based 
upon  a  cubic  yard  of  plastic  concrete;  hence,  the  quantities  were  scaled 
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down  for  use  in  the  laboratory.  The  estimated  quantities  were  then  mixed 
and  the  properties  of  the  fresh  (slump  and  workability)  and  hardened 
concrete  (strength)  observed.  The  first  trial  mix  did  not  yield  the  de¬ 
sired  properties  of  slump,  workability,  nor  strength.  Adjustments  were 
made  in  the  water-cement  ratio  and  in  the  relative  proportions  of  the 
sand  and  gravel.  The  water-cement  ratio  adjusted  the  strength  and  the 
ratio  of  sand  to  gravel  affected  the  slump  and  workability. 

After  five  or  six  trial  mixes  and  adjustments,  a  concrete  mix  was 
obtained  that  produced  the  desired  slump,  workability,  and  strength.  The 
relative  proportions  by  weight  of  the  final  concrete  mix  before  free  mois¬ 
ture  or  absorption  adjustments  were: 

Water-cement  ratio  0.56 

Cement-sand-gravel  1  -  2.36  -  2.45 

This  concrete  mix  had  a  rather  large  ratio  of  sand  to  gravel.  The 
gravel  was  angular  shaped  and  therefore  required  a  greater  amount  of  sand 
to  obtain  the  desired  workability. 

2.4  MIXING  PROCEDURE 

The  water,  cement,  sand,  and  gravel  were  mixed  in  a  power-driven 
revolving  drum  mixer  that  had  a  capacity  of  approximately  two  cubic  feet. 
Prior  to  mixing,  a  mortar  composed  of  equal  parts  of  cement  and  sand  was 
placed  in  the  mixer.  The  mixer  was  then  run  for  a  few  minutes  until  the 
mortar  had  covered  the  entire  wall  of  the  mixer.  The  excess  was  then 
discarded.  The  mortar  adhering  to  the  mixer  was  intended  to  compensate 
for  loss  of  mortar  from  the  batch. 
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The  ingredients  were  weighed  to  the  nearest  one-tenth  of  a  pound. 

The  gravel,  sand  and  cement  were  placed  in  the  mixer  in  the  order  men¬ 
tioned.  The  water  was  slowly  added  while  the  mixer  was  running.  The 
mixing  continued  for  approximately  five  minutes  after  all  the  water  had 
been  added. 

A  slump  test  was  made  immediately  after  mixing.  The  slump  tests 
were  conducted  according  to  ASTM  standards.  The  results  of  all  tests 
are  shown  in  Table  1.  The  desired  slump  was  3-4  inches.  This  range  was 
not  always  obtained  as  noted  from  Table  1.  The  slump  was  quite  sensitive 
to  slight  changes  in  relative  percentages  of  different  particle  sizes  and 
to  humidity  and  temperature  in  the  laboratory. 

After  the  slump  test  had  been  completed,  the  concrete  was  placed  in 
the  cylinder  and  cube  molds.  The  cylinders  were  the  standard  6  x  12  inch 
test  cylinders  and  were  used  for  quality  control.  The  cube  molds  were 
made  of  aluminum  plate  1/2  inch  thick.  A  mold  is  shown  in  Figure  2. 

There  were  six  cube  molds  per  unit.  The  nominal  inside  mold  dimensions 
were  3  inches  on  each  side.  The  concrete  cylinders  were  cast  according 
to  ASTM  standards. 

Consolidation  of  the  concrete  in  the  cube  molds  was  accomplished  by 
vibrating  the  molds  while  the  concrete  was  being  placed  in  the  molds. 

The  concrete  was  placed  in  the  molds  in  approximately  three  equal  layers. 
Vibration  was  continued  only  long  enough  to  obtain  a  relatively  smooth 

i 

surface  of  the  concrete.  Different  methods  of  consolidating  the  concrete 
in  the  cube  molds  had  been  attempted.  Tamping  with  different  size  rods 
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TABLE  1 


SLUMPS  FOR  CONCRETE  MIXES 


Date 

Slump  (in) 

1-28-71 

5.75 

2-1-71 

4.00 

2-2-71 

2.75 

2-3-71 

3.24 

2-4-71 

3.24 

2-8-71 

2.50 

2-9-71 

3.25 

2-11-71 

3.00 

2-16-71 

3.25 

2-24-71 

3.00 

2-25-71 

3.50 

3-2-71 

3.75 

3-8-71 

3.75 

3-10-71 

4.00 

3-11-71 

4.00 

3-15-71 

5.00 

3-17-71 

5.00 

3-22-71 

5.50 

3-29-71 

3.75 

3-31-71 

3.50 

4-7-71 

3.50-3.00 

4-12-71 

4.00-2.75 

4-14-71 

3.75-6.00 

4-19-71 

3.50-4.50 

4-21-71 

3.75-4.00 

4-26-71 

2.75-3.25 

4-28-71 

3.50-3.50 

4-29-71 

3.75-3.75 

5-3-71 

4.00-5.75 

5-5-71 

4.50-5.00 

5-6-71 

4.50-4.50 

5-10-71 

3.50-3.75 

5-12-71 

4.50-4.50 

5-19-71 

4.00-3.75 

14 


Figure  2.  Cube  Molds 
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was  attempted.  Also  tapping  the  mold  and  tamping  simultaneously  was 
attempted.  The  effectiveness  of  the  different  methods  of  consolidation 
was  determined  by  sawing  the  cubes  into  halves  and  visually  noting  the 
number  and  size  of  voids  in  the  hardened  concrete.  The  method  described 
above  yielded  the  smallest  number  and  smallest  sized  voids. 

The  cylinder  molds  were  coated  with  oil  and  the  cube  molds  were 
coated  with  wheel  bearing  grease  prior  to  pouring  of  the  concrete.  This 
was  to  prevent  the  concrete  from  sticking  to  the  walls  of  the  molds.  A 
thin  layer  of  the  grease  was  used. 

The  specimens  were  allowed  to  set  in  the  molds  for  about  24  hours. 
After  24  hours,  the  specimens  were  removed  from  the  molds  and  placed  in 
a  curing  room  for  28  days.  The  curing  room  conditions  were  100  percent 
humidity  and  approximately  75  degrees  F. 

The  cylinders  were  tested  for  compressive  strength  at  7  days  and  at 
28  days.  The  results  of  these  tests  are  shown  in  Table  2.  A  noticeable 
increase  in  cylinder  strength  occurred  in  the  cylinders  cast  on  March  8 
and  thereafter.  This  increase  in  strength  was  thought  to  have  been  due 
to  the  use  of  a  new  cement.  The  new  cement  was  of  the  same  type  as  the 
cement  used  previously. 

\  * 

A  typical  strength  curve  in  which  the  strength  is  plotted  versus 
curing  time  is  shown  in  Figure  3.  This  curve  was  obtained  by  casting 
14  cylinders  using  the  same  mix  and  testing  3  cylinders  at  7,  14,  and 
28  days.  Two  cylinders  were  tested  at  35  days  of  curing.  All  cylinders 
including  the  ones  tested  at  35  days  were  moist  cured  until  the  day  of 
testing. 
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TABLE  2 


CYLINDER  STRENGTHS 


7 -day 

strength 

28-day  strength 

Date  Mixed 

lbs. 

Psi 

lbs. 

Psi 

2-2-71 

90,500 

3200 

118,000 

4173 

2-3-71 

83,500 

2953 

122,000 

4314 

2-4-71 

94,000 

3324 

122,500 

4332 

2-8-71 

98,500 

3483 

128,000 

4527 

2-9-71 

91,500 

3236 

125,000 

4420 

2-11-71 

78,000 

2758 

113,000  (35  day) 

3996  (35  day) 

2-16-71 

90,000 

3183 

122,500 

4332 

2-24-71 

80,000 

2829 

118,000 

4173 

2-25-71 

80,000 

2829 

116,000 

4102 

3-2-71 

73,500 

2599 

114,500 

4049 

3-8-71 

92,500 

3271 

136,500 

4827 

3-10-71 

95,000 

3359 

136,000 

4810 

3-11-71 

87,500 

3094 

136,000 

4810 

3-15-71 

91,000 

3218 

132,500 

4686 

3-17-71 

96,000 

3395 

140,500 

4969 

3-22-71 

90,000 

3183 

134,000 

4739 

3-29-71 

95,500 

3377 

130,000 

4597 

3-31-71 

89,000 

3147 

131,000 

4633 

4-7-71 

100,500 

3554 

137,000 

4845 

4-12-71 

94,000 

3324 

140,000 

4951 

4-14-71 

101,000 

3572 

143,000 

5057 

4-19-71 

96,000 

3392 

129,500 

4580 

4-21-71 

89,000 

3147 

129,500 

4580 

4-26-71 

77,000 

2723 

118,500 

4191 

4-28-71 

75,500 

2670 

117,000 

4138 

4-29-71 

88,500 

3130 

144,500 

5110 

5-3-71 

90,500 

3200 

141,500 

5004 

5-5-71 

89,000 

3147 

114,500 

4058 

5-6-71 

93,500 

3306 

143,000 

5057 

5-10-71 

90,500 

3200 

139,000 

4916 

5-12-71 

92,500 

3271 

128,500 

4544 

5-19-71 

88,500 

3130 

120,000 

4244 

5-24-71 

104,000 

3678 

142,000 

5013 

17 
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Figure  3.  Strength  Curve  for  Concrete 


SECTION  III 


DESCRIPTION  OF  EQUIPMENT 


3.1  TESTING  MACHINE 

The  triaxial  testing  machine  was  designed  and  constructed  as  a  part 
of  the  research  program.  The  testing  machine  design  was  based  upon  the 
following  criteria: 

a.  Loads  could  be  applied  in  three  orthogonal  directions. 

b.  The  loading  in  each  direction  would  be  independently  controlled. 

c.  The  loading  capacity  of  the  testing  machine  would  be  such  that 
a  maximum  stress  of  30,000  psi  could  be  applied  triaxially  to 
a  three-inch  cube. 

d.  The  measurement  of  the  applied  load  and  the  accompanying  deforma¬ 
tions  would  be  possible. 

e.  The  frame  would  deform  symmetrically  under  an  applied  load. 

f.  The  loads  could  be  applied  in  compression  or  tension. 

3.2  TESTING  FRAME 

The  frame  of  the  testing  machine  was  dimensioned  such  that  there 
was  sufficient  space  for  the  hardware  such  as  the  jacks,  compression 
platens,  load  cells,  and  extensometers  to  be  installed. 

The  testing  frame  was  designed  to  be  flexible  and  to  deform  sym¬ 
metrically  such  that  the  test  specimen  would  remain  in  a  fixed  position 
as  the  loads  were  applied.  Upon  checking  the  performance  of  the  machine, 
it  was  noted  that  a  slight  shift  of  the  test  specimen  occurred  as  the 
loads  were  applied.  The  shift  was  in  a  direction  away  from  the  active 


19 


ja<ks.  A  shift  in  the  Lest  cube  would  not  affect  the  uniaxial  tests. 

It  would  introduce  friction  forces  between  the  Lest  cube  and  the 
loading  platens  in  biaxial  and  triaxial  tests.  Friction  reducing 
methods  were  utilized  in  the  testing  program;  hence,  the  magnitude 
of  the  induced  friction  forces  was  probably  small. 

The  main  members  of  the  testing  frame  consisted  of  wide-flanged 
steel  sections.  The  sections  were  selected  according  to  AISC  specifica¬ 
tions  for  A36  steel.  A  load  factor  of  two  was  also  applied  in  the 
structural  design. 

All  connections  on  the  testing  trames  were  welded.  The  spans  were 
short  and  the  possible  loads  high;  hence,  a  large  number  of  shear  stif¬ 
fener  plates  were  required.  Web  stiffners  were  also  used  beneath  the 
jacks  and  the  bearing  plates  to  prevent  buckling  of  the  webs. 

The  testing  machine  was  composed  of  two  frames,  one  being  suspended 
horizontally  within  the  other  as  shown  in  Figure  4.  The  horizontal 
frame  was  suspended  from  the  vertical  frame  by  two  cable  hoists.  Two 
jacks  were  attached  to  the  horizontal  frame  and  one  to  the  vertical 
frame.  The  horizontal  frame  was  adjustable  in  both  the  horizontal  and 
vertical  directions. 

The  loads  were  applied  by  means  of  manually  operated  hydraulic  jacks 
with  capacities  of  150  tons.  The  jacks  were  equipped  with  spherical 
heads  which  were  self  aligning.  The  jacks  were  also  double-acting. 

This  was  necessary  in  order  that  tensile  loads  could  be  applied. 

The  compression  platens  were  machined  from  high  strength  steel  (75 
ksi).  The  bearing  surface  dimensions  were  2.75  inches  by  2.75  inches. 
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The  area  of  the  platens  were  necessarily  smaller  than  the  tesc  specimens 
to  prevent  the  platens  from  contacting  against  each  other  during  a  test. 
Tne  compression  platens  had  a  slot  milled  in  them  to  accomodate  the 
extensometer  arms.  In  addition,  a  hole  was  drilled  into  the  center  of 
the  bearing  surface  to  accomodate  a  flange  pin  as  shown  in  Figures  5  and 
6.  The  extensometer  arm  was  placed  in  the  slot  and  on  the  one  end  of  the 
pin.  The  purpose  of  the  pin  was  to  remove  deformations  arising  from  the 
use  of  friction  pads  such  that  the  deformation  detected  by  the  extenso- 
meters  was  the  deformation  at  the  center  of  the  test  cube  only.  Whenever 
friction  pads  were  used,  polythelene  washers  of  the  same  thickness  as 
the  friction  reducing  pads  were  placed  behind  the  flange  of  the  pins  as 
shown  in  Figure  6.  The  friction  reducing  pads  had  a  hole,  which  was 
slightly  larger  in  diameter  than  the  diameter  of  the  pin  flange,  punched 
in  their  centers.  This  allowed  the  pin  flange  to  directly  contact  the 
surface  of  the  test  cube;  hence,  the  only  deformation  detected  was  that 
of  the  center  of  the  cubes. 

A  pair  of  tension  heads  are  shown  in  Figure  7.  The  tension  heads 
are  connected  to  the  test  specimen  with  an  epoxy  glue  and  to  the  load 
cells  by  means  of  ball  joints.  All  of  the  tension  heads  had  a  1/4  inch 
hole  drilled  at  their  centers  to  accomodate  steel  pins.  The  extenso¬ 
meter  arms  were  then  connected  to  these  Tins.  The  contact  area  of  the 
tension  heads  was  three  inches  square. 

3.3  LOAD  CELLS 

The  load  cells  used  in  the  testing  program  were  constructed  and 
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Figure  5.  Compression  Platen 
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Figure  7.  Tension  Heads 
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calibrated  at  NMSU.  The  cells  consisted  of  columns  five  inches  in 
length,  made  of  steel  or  aluminum.  Four  electrical  resistance  strain 
gages  were  mounted  at  mid-height  of  the  columns.  The  strain  gages  were 
arranged  and  wired  in  a  manner  such  that  only  axial  strains  were 
detected. 

Three  different  types  of  load  cells  were  used  for  the  entire  test¬ 
ing  program.  Two  types  were  used  in  compression  testing  and  a  third 
type  was  used  in  tension  testing.  The  three  types  of  load  cells  are 
shown  in  Figure  8.  One  type  of  compression  load  cells  had  circular 
cross-sections  and  was  made  of  aluminum.  Their  capacity  was  120,000 
pounds.  The  second  type  of  compression  load  cells  had  three-inch  square 
cross-sections  and  were  made  of  high-strength  steel.  Their  capacity 
was  300,000  pounds.  The  aluminum  load  cells  were  used  in  the  biaxial 
tests  and  the  steel  cells  in  the  triaxial  tests.  Two  types  of  compres¬ 
sion  load  cells  were  used  to  obtain  better  accuracy  at  the  lower  load 
levels. 

The  tension  load  cells  had  circular  cross-sections  and  were  made 
of  steel.  Their  capacity  was  controlled  by  the  one-half  inch  attaching 
bolts  whose  capacity  was  approximately  5,000  pounds. 

The  load  cells  were  calibrated  with  the  aid  of  a  universal  testing 
machine  located  in  the  Materials  Testing  Laboratory  at  NMSU.  The  maximum 
capacity  available  was  200,000  pounds.  This  exceeded  the  capacity  of 
the  aluminum  compression  load  cells  and  the  tension  load  cells,  but  was 
less  than  the  capacity  of  the  steel  compression  load  cells.  The  300,000 
pound  cells  were,  therefore,  calibrated  by  loading  up  to  200,000  pounds. 
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Figure  8.  Load  Cells 
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It  was  assumed  that  the  load-strain  relationship  was  valid  up  to  300 >000 
pounds . 

The  same  calibration  procedure  was  followed  for  each  type  of  load 
cell.  The  cells  were  loaded  incremently  with  load  and  strain  readings 
taken  after  each  increment.  About  ten  readings  were  taken  between  no 
load  and  maximum  load.  The  cells  were  then  unloaded  incremently  and 
readings  taken  again  after  each  unloading  increment.  This  procedure 
was  repeated  until  readings  were  taken  for  five  loading  and  unloading 
cycles.  All  of  the  load  and  strain  readings  were  then  taken  and  a 
straight  line  relationship  was  fitted  using  the  method  of  least  squares. 
The  calibration  constants  for  the  load  cells  are  given  below: 


LOAD  CELLS 

Units 

CELL  NUMBERS 
(lbs/micro  in.  per 

inch) 

No.  1 

No.  2 

No.  3 

Aluminum 

15.790 

15.295 

15.920 

Steel 

(Compression) 

89.430 

88.060 

89.510 

Steel 

(Tension) 

6.270 

6.803 

7.076 

3.4  EXTENSOMETERS 

The  extensometers  consist  of  an  aluminum  bar  attached  to  two  steel 
clamps  shaped  as  shown  in  Figure  9.  The  elbows  of  the  steel  clamps 
are  necessary  to  hold  the  extensometer  away  from  the  loading  platens  on 
the  other  loading  axes.  An  extensometer  in  position  for  testing  is 
shown  in  Figure  10.  Four  electrical  resistance  strain  gages  were 
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Figure  9.  Extensometer 
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View  of  Extensometers  in  Place 


attached  to  the  aluminum  bar.  They  were  arranged  and  wired  in  a  manner 
such  that  only  flexural  strains  were  detected. 

The  extensometers  were  calibrated  in  the  testing  frame  and  in  the 
same  position  as  used  in  testing.  This  was  found  to  be  necessary  when 
it  was  discovered  that  a  different  calibration  constant  resulted  if  a 
given  extensometer  was  used  in  different  positions  and  on  different  axes. 

Calibration  of  the  extensometers  was  achieved  with  the  aid  of  a 
hexagonal  shaped  aluminum  bar.  The  bar  was  three  inches  in  length  and 
had  two  electrical  resistance  strain  gages  mounted  at  its  mid-height. 

The  strain  gages  were  wired  to  detect  axial  strains  only.  The  aluminum 
bar  was  placed  into  the  testing  machine  between  two  compression  platens 
of  a  particular  axis.  The  strain  gages  were  wired  into  an  X-Y  recorder. 
The  extensometer  strain  gages  were  wired  into  the  other  axis  of  the  X-Y 
recorder.  A  load  was  then  applied  to  the  aluminum  bar.  The  result 
indicated  on  the  X-Y  recorder  graph  is  the  strain  indicated  by  the 
aluminum  bar  versus  the  strain  output  of  the  extensometer.  The  aluminum 
bar  was  three  inches  in  length,  the  same  as  the  nominal  length  of  the 
concrete  test  specimens;  hence,  it  was  assumed  that  the  strains  indicated 
by  the  aluminum  bar  would  be  comparable  to  the  strains  in  the  test  speci¬ 
mens.  The  calibration  constant  was  obtained  by  taking  the  average 
slope  of  the  plot  obtained  on  the  X-Y  recorder.  A  typical  record  of 
this  type  is  shown  in  Figure  11.  The  load  was  applied  and  released  a 
number  of  times  to  get  an  indication  of  the  extensometers  ability  to 
duplicate  strains.  This  is  indicated  by  the  width  of  the  band.  The 
extensometers  were  removed  and  replaced  during  calibration  to  determine 
if  this  changed  the  constant.  Removing  and  replacing  the  extensometers 
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did  not  appear  to  change  the  calibration  constant  as  did  a  change  of 
axis  and  position. 

The  extensometers  were  recalibrated  whenever  any  hardware  such  as 
load  cells  or  platens  were  changed.  A  change  in  hardware  did  usually 
result  in  a  slight  change  in  the  calibration  constant. 

3.5  RECORDING  EQUIPMENT 

The  load  and  deformation  were  monitored  for  each  of  the  three  axes. 
The  deformation  was  measured  at  the  center  of  the  concrete  cube.  This 
information  was  recorded  with  the  use  of  X-Y  recorders.  To  accomplish 
this,  amplifiers  were  used  to  obtain  the  desired  accuracy  and  flexibility 
in  changing  scales  for  the  different  tests.  The  loads  and  deformations 
were  scaled  such  that  the  recorded  information  was  the  stress-strain 
curve  for  the  given  axis.  The  loads  were  detected  using  load  cells  and 
the  deformations  were  detected  using  special  extensometers. 

The  X-Y  recorders  used  were  manufactured  by  Instron,  Mosley,  and 
EAI.  The  amplifiers  were  manufactured  by  Instron,  Hallmark,  and  Newport. 
The  arrangement  of  the  instrumentation  is  shown  in  Figure  12. 

Both  the  vertical  and  horizontal  axes  of  the  recorders  were  cali¬ 
brated  using  a  precision  fixed  resistor  of  known  value.  The  resistor 
most  often  used  produced  a  strain  of  518  p  in/in  when  placed  across  one 
arm  of  a  four-arm  bridge. 

3.6  COMMENTS  REGARDING  THE  TESTING  MACHINE 

An  attempt  was  made  to  evaluate  the  performance  of  the  testing 
machine.  Six  electrical  resistance  strain  gages  were  mounted  (three 
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Figure  12.  Test  Equipment 
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parallel  to  each  other  on  two  opposite  sides)  on  a  three- inch  aluminum 
cube.  The  aluminum  cube  was  then  placed  in  the  testing  machine  and 
loaded  uniaxially.  The  load  was  applied  parallel  and  perpendicular  to 
the  directions  of  the  strain  gages.  A  calibrated  load  cell  was  used  to 
determine  the  magnitude  of  the  applied  load  and  a  calibrated  extensometer 
was  used  to  determine  the  deformation  of  the  center  of  the  cube.  One 
purpose  of  the  test  was  to  determine  if  the  extensometer  could  be  used 
to  determine  strain  in  the  cube. 

Strain  readings  were  taken  to  determine  the  distribution  of  strains 
throughout  the  cube.  A  uniform  strain  distribution  proved  difficult  to 
obtain.  A  slight  eccentricity  would  cause  rather  large  differences  in 
strain  from  side  to  side.  Finally,  a  transit  was  used  to  align  the 
jack,  compression  platens,  cube  and  load  cell.  Even  when  aligned  with 
the  transit,  rather  large  variations  in  strain  were  obtained.  The  strain 
readings  were  also  compared  to  the  readings  obtained  with  the  extenso¬ 
meter.  It  was  noted  that  the  average  of  the  strain  readings  taken  from 
the  gages  mounted  on  the  aluminum  cube  compared  quite  favorably  with 
the  values  obtained  with  the  extensometer.  It  was,  therefore,  concluded 
that  the  extensometers  could  be  used  in  determining  strains  in  concrete 
cubes.  This  was  also  verified  after  uniaxial  test  results  on  concrete 
became  available.  The  stress-strain  curves  obtained  from  the  tests  were 
comparable  to  those  published  in  various  journals. 
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SECTION  IV 


TESTING  PROCEDURES 


4.1  INTRODUCTION 

The  procedures  utilized  during  the  testing,  the  problems  encountered 
during  testing,  and  the  problem  solutions  or  compromises  are  reviewed  in 
this  section. 


4.2  TEST  SPECIMENS 

4.2.1  Compression  Tests 

The  three  inch  cubical  concrete  test  specimens  were  cast  in  aluminum 
molds  as  described  in  Section  II.  Five  of  the  six  sides  of  the  cubes 


were  formed  by  the  sides  of  the  molds.  The  sixth  side  was  trowelled 
after  the  concrete  had  been  consolidated  in  the  molds;  hence,  it  was 
not  nearly  so  smooth  nor  as  plane  as  the  other  sides  of  the  cubes.  The 
trowelled  sides  were  later  trimmed  with  a  diamond  saw  blade  and  then 
sanded  on  a  disc  sander  to  make  it  smoother.  This  surface  will  be  re¬ 
ferred  to  as  the  shaped  surface. 

Various  surface  treatment  methods  were  considered.  They  are  listed 
below.  The  maximum  difference  in  cube  dimensions  from  point  to  point 
for  the  different  methods  considered  are  also  indicated  below. 


Surface  Treatment 
No  treatment 
Lightly  sanded 
Well  sanded 


Dimensional  Variation  (inches) 


Formed  Surfaces 

Shaped  Surface 

0.008 

0.020 

0.011 

0.013 

0.015 

0.018 
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The  above  values  do  not  give  an  indication  of  flatness  or  plane¬ 
ness  of  the  surfaces.  The  flatness  was  more  critical  than  the  difference 
in  dimensions  since  adjustable  heads  were  used  on  the  jacks. 

Unconfined  compression  tests  (uniaxial)  were  conducted  to  determine 
any  effects  due  to  the  different  surface  treatments.  The  results  are 
indicated  below.  The  listed  strength  values  were  obtained  by  averaging 
results  from  three  tests. 

Surface  Treatment  Avr.  Uniaxial  Strength 

(a)  Formed  surfaces  lightly  sanded, 

shaped  surface  trimmed  and 
sanded,  loaded  on  shaped  sur¬ 
face.  4270  psi 

(b)  Formed  surfaces  lightly  sanded, 

shaped  surface  trimmed  and 
sanded,  loaded  on  sanded  sur¬ 
face.  3850  psi 

(c)  All  surfaces  trimmed  and  sanded, 

loaded  on  shaped  surface.  3780  psi 

(d)  All  surfaces  trimmed  and  sanded, 

loaded  on  sanded  surfaces  3840  psi 

From  the  above  results  it  was  concluded  that  the  surfaces  could  be 
lightly  sanded  and  loaded  on  the  sanded  surfaces.  This  treatment  was 
the  simplest.  Treatments  (b) ,  (c),  and  (d)  were  not  considered  to  be 
significantly  different.  At  that  point  in  time  it  had  not  been  recognized 
that  the  difference  in  strengths  obtained  from  loading  on  the  sanded  and 
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shaped  surfaces  was  not  necessarily  due  to  surface  treatment,  but  due 
to  the  anisotropic  cubes  resulting  from  the  casting  operations. 

The  uniaxial  and  biaxial  tests  were  conducted  on  cubes  that  were 
lightly  sanded  on  the  formed  sides  and  trimmed  on  the  trowelled  sides. 
Loads  for  these  tests  were  applied  to  the  sanded  surfaces  and  not  to 
the  trowelled  surfaces. 

Later  for  some  of  the  triaxial  tests,  the  trowelled  sides  were 
trimmed  with  the  diamond  blade  and  then  placed  into  a  surface  grinding 
machine.  The  finished  surface  was  then  quite  smooth  and  plane.  The 
maximum  dimensional  variation  was  usually  less  than  0.003  inch  whenever 
the  surface  grinder  was  used. 

Np 

Uniaxial  tests  on  cubes  whose  trowelled  side  was  finished  in  the 
surface  grinder  indicated  that  the  difference  in  strengths  from  test  to 
test  was  smaller  than  for  the  cases  where  the  trowelled  surface  was  not 
planed. 


4.2.2  Tension  Tests 

No  special  surface  treatment  was  required  in  the  case  of  cubes  sub¬ 
jected  to  tensile  loads.  Tension  heads  were  attached  to  the  cubes  by 
means  of  an  epoxy  resin  layer  which  was  approximately  1/8  inch  in  thick¬ 
ness.  The  tensile  strength  of  the  epoxy  resin  was  much  larger  than  the 
tensile  strength  of  the  concrete.  The  modulus  of  elasticity  of  most  un¬ 
filled  epoxies  (no  fillers)  ranges  from  50,000  psi  to  500,000  psi  as 
compared  to  a  modulus  of  elasticity  of  3-6  million  psi  for  concrete. 

The  low  modulus  of  elasticity  of  the  epoxy  coupled  with  the  1/8  inch 
lay  would  provide  little  restraint  at  the  tension  heads. 
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4.3  UNIAXIAL  COMPRESSIVE  TESTS 


Uniaxial  test  results  were  frequently  used  to  normalize  results 
obtained  from  biaxial  and  triaxial  tests;  hence,  a  reliable  uniaxial 
strength  value  was  desired.  Only  18  cube  molds  were  available;  hence, 
a  rather  large  number  of  batches  of  concrete  was  required  to  obtain  a 
sufficient  number  of  test  cubes.  There  was  batch  to  batch  variation  in 
ultimate  strengths.  This  made  it  necessary  to  perform  uniaxial  tests 
for  each  batch. 

The  uniaxial  tests  were  conducted  in  the  triaxial  testing  machine. 
The  compression  platens  on  the  testing  machine  were  smaller  in  area  than 
the  area  of  the  cubes.  In  addition  the  cubes  were  beveled  along  their 
edges.  This  also  reduced  the  cross-sectional  area.  The  length  to  width 
ratio  for  cubes  is  unity;  hence,  friction  effects  between  the  compres¬ 
sion  platens  and  cubes  would  greatly  affect  the  uniaxial  test  results. 


4.3.1  Friction 

Several  uniaxial  tests  were  performed  using  friction  reducing  pads 
and  several  tests  without  friction  reducers.  The  friction  reducing  pads 
consisted  of  two  polyethelene  sheets  four  mils  thick  with  grease  between 
the  sheets  as  suggested  by  Mills  (23).  Other  friction  reducing  methods 
were  tried;  however,  the  method  mentioned  above  yielded  thu  lowest  uni¬ 
axial  strengths.  The  results  of  these  tests  are  shown  beJow. 


No  Friction  Reducer 
5760  psi 
4990  psi 
5080  psi 
Avg.  5280  psi 


Friction  Reducei 
3530  psi 
3830  psi 
3660  psi 
Avg.  3670  psi 
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The  results  shown  on  the  previous  page  clearly  indicate  that 
friction  effects  are  considerable  in  uniaxial  tests  of  cubes.  As  a 
result  the  friction  reducing  pads  consisting  of  the  two  polyethelene 
(a  thermoplastic  synthetic  resin)  sheets  and  grease  were  used  in  all 
tests. 


4.3.2  Beveling 

The  concrete  test  cubes  were  beveled  along  their  edges.  This  was 
performed  to  reduce  the  contact  area  of  the  cubes  to  approximately  the 
same  area  as  the  compression  platens.  Beveling  also  reduced  the  cross- 
sectional  area  of  the  cubes.  Uniaxial  tests  were  conducted  to  deter¬ 
mine  the  effects  due  to  beveling.  Beveled  and  unbeveled  cubes  were 
tested  in  a  universal  testing  machine.  Friction  reducing  pads  were  used 
in  these  tests.  The  results  of  these  tests  were  compared  to  the  re¬ 
sults  obtained  by  stacking  four  cubes  and  testing  uniaxially  in  a  uni¬ 
versal  testing  machine.  This  comparison  was  considered  of  interest 
since  in  the  stacked  cube  tests,  friction  effects  are  essentially 
eliminated  as  a  result  of  the  larger  height  to  width  ratio  and  not  by 
the  use  of  friction  reducing  pads.  The  results  of  the  tests  are  shown 
below.  In  all  cases  the  stresses  were  computed  on  the  basis  of  a 
three-inch  square  cross  section. 


Cubes 


Batch  No.  1 


Batch  No.  2 


Stacking  Four  Cubes 
Cubes  Not  Beveled 
Cubes  Beveled 


4720  psi 
4600  psi 
4040  psi 


4390  psi 
4570  psi 
4010  psi 
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The  percentage  difference  between  .he  beveled  and  unbeveled  cubes 
for  the  two  batches  was  -12.2  and  -12.4,  respectively.  The  difference 
between  results  of  the  stacked  cubes  and  the  unbeveled  cubes  was  not 
considered  to  be  significant.  The  differences  were  no  more  than  the 
test  to  test  variation. 

4.3.3  Testing 

Whenever  conducting  uniaxial  tests  in  the  triaxial  testing  machine, 
care  was  taken  to  ensure  that  no  lateral  confinement  would  occur  until 
at  least  the  maximum  stress  had  been  obtained.  The  slightest  lateral 
pressure  would  greatly  affect  the  test  results.  The  compressive  platens 
were  placed  adjacent  to  the  unconfined  sides  of  the  test  cube.  This  was 
necessary  whenever  lateral  strains  were  measured.  The  load  cells  were 
activated  in  the  directions  of  the  unconfined  sides.  Any  pressure  that 
developed  could  therefore  be  detected.  A  spacer  was  used  to  prevent 
contact  initially  with  the  platens  located  near  the  unconfined  sides  to 
prevent  lateral  pressure  from  developing  during  a  test. 

4.4  UNIAXIAL  TENSION  TESTS 

There  were  no  significant  problems  involved  with  the  uniaxial  ten¬ 
sion  tests.  The  tension  heads  which  were  glued  to  the  cubes  with  an 
epoxy  resin  were  connected  to  the  load  cell  and  jack  head  by  means  of 
ball  joints.  The  ball  joints  allowed  the  cube  to  align  itself  as  the 
tensile  load  was  applied.  Special  hardware  was  required  to  protect  the 
extensometers  whenever  separation  of  the  cubes  occurred.  Strains  were 
small  in  the  uniaxial  tension  tests;  hence,  the  accuracy  of  the  recorded 
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strains  may  be  questionable.  The  extensometers  were  designed  to  detect 
large  strains. 

During  the  uniaxial  tension  test,  it  was  discovered  that  different 
strength  values  could  be  obtained  depending  upon  the  direction  the  ten¬ 
sile  load  was  applied.  Higher  strength  values  were  obtained  whenever 
the  load  was  applied  to  the  formed  sides  of  the  cubes.  The  results  of 
four  uniaxial  tension  tests  are  shown  below. 

Load  Applied  to  Shaped  Load  Applied  to  Formed 

_ Surface _ Surfaces _ 

237  psi  364  psi 

287  psi  322  psi 

Avg.  262  psi  Avg.  343  psi 

The  cubes  were  all  from  the  same  batch.  The  percentage  differ¬ 
ences  between  the  average  values  is  31.  This  indicates  that  the  concrete 
cubes  were  not  isotropic.  The  properties  are  essentially  the  same  in 
the  directions  of  the  cubes  formed  by  the  sides  of  the  molds,  but  are 
somewhat  different  in  the  direction  of  casting. 

4. 5  BIAXIAL  TESTS 

4.5.1  Compressive  Tests 

Biaxial  tests  were  conducted  to  determine  the  effect  of  using  fric¬ 
tion  pads.  As  in  the  case  of  the  uniaxial  tests,  the  test  results  with¬ 
out  the  use  of  the  friction  reducing  pads  indicated  higher  strengths. 

In  addition,  the  tests  without  friction  pads  indicated  a  more  ductile 
material  than  did  the  tests  where  friction  reducing  pads  were  used; 
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therefore,  all  biaxial  tests  were  conducted  using  friction  reducing 
pad  s 

Biaxial  test  results  were  also  greatly  affected  by  the  presence  of 
a  slight  lateral  pressure  on  the  un^onfined  axis;  therefore,  a  small 
space  was  left  between  the  cube  and  the  compression  platens.  The  load 
cell  in  that  direction  was  activated;  hence,  the  presence  of  any  lateral 
pressure  was  easily  detected. 

The  biaxial  tests  were  conducted  using  the  jacks  on  the  horizontal 
frame  of  the  triaxial  testing  machine.  The  test  cubes  were  always 
placed  into  the  testing  machine  in  the  same  manner,  that  is,  the  orienta¬ 
tion  was  the  same  as  in  the  casting  operation.  The  shaped  surface  faced 
upward.  No  load  was  applied  to  the  shaped  surface  in  the  biaxial  test 
series  until  late  in  the  testing  program.  In  one  test  the  major  load 
was  applied  to  the  shaped  surface.  The  strength  values  were  larger  than 
had  been  previously  obtained.  A  test  series  was  then  initiated  to  de¬ 
termine  if  the  orientation  of  loading  affected  the  test  results.  The 
test  results  are  shown  below.  The  values  shown  are  averages  of  two 
tests.  Equal  loads  were  applied  in  the  two  directions. 

(a)  Load  Applied  to  Formed  Surfaces  4040  psi 

(b)  Loads  Applied  to  Two  Formed  Surfaces 

and  the  Shaped  Surface  5300  psi 

There  is  a  rather  large  difference  between  the  strengths  resulting 
from  applying  the  loads  to  the  formed  surfaces  or  to  the  shaped  surir>ce. 
Fortunately  all  of  the  cubes  were  loaded  the  same  way  in  the  biaxial 
tests. 
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The  first  few  biaxial  tests  were  performed  by  applying  the  minor 
load  and  then  applying  the  major  load  until  failure  occurred.  This 
loading  procedure  was  abandoned  after  only  a  few  tests.  It  was  very 
difficult  to  control  the  minor  load  at  a  constant  value  during  the  time 
the  major  load  was  being  applied. 

Proportional  loading  (load  in  each  direction  applied  simultaneously 
and  in  a  fixed  proportion)  was  used  for  the  majority  of  the  biaxial 
tests.  The  loads  were  easier  to  control  whenever  proportional  loading 
was  used.  In  addition,  the  results  from  a  proportional  loading  test 
were  easier  to  interpret. 

Incremental  loading  was  used  on  a  few  biaxial  tests  for  the  pur¬ 
pose  of  studying  the  interaction  between  the  two  loaded  axes. 

The  time  required  to  fail  a  cube  subjected  to  proportion  biaxial 
loads  was  less  than  two  minutes. 

Once  a  maximum  stress  had  been  reached,  it  was  impossible  to 
maintain  the  same  proportions  between  the  applied  load;  hence,  the 
loading  proportions  listed  on  any  of  the  stress-strain  records  was 
valid  only  up  to  the  maximum  stress  on  the  major  axis. 

Some  difficulty  was  encountered  in  conducting  tests  in  which  the 
ratio  of  the  major  axis  load  to  the  minor  axis  load  was  greater  than  10. 
The  application  of  the  major  load  would  produce  a  load  in  the  minor  axis 
that  was  greater  than  the  intended  value.  This  required  a  release  of 
pressure  in  the  jack  applying  the  minor  load.  Releasing  the  pressure 
and  maintaining  a  constant  load  level  was  difficult. 


4.5.2  Tension-Tension  and  Compressive-Tension  Tcs Is 
The  greatest  difficulty  encountered  in  the  biaxial  tension  and  com¬ 
pressive-tension  tests  was  the  elimination  of  bending  stresses  induced 
durin  -  testing.  This  problem  was  never  resolved.  An  initially  per  ¬ 
fectly  aligned  test  specimen  would  deform  as  the  loads  \.ere  applied  such 
that  bending  stresses  were  induced  into  the  test  specimen.  In  the  case 
of  compressive-tension  tests,  the  larger  the  compressive  force  in  rela¬ 
tion  to  the  tensile  force,  the  greater  the  induced  bending  stresses. 

It  is  believed  that  the  indicated  maximum  stresses  in  this  loading 
region  are  less  than  the  actual  concrete  strength.  The  error  intro¬ 
duced  by  the  addition  of  bending  is  difficult  to  estimate.  The  magni¬ 
tude  of  the  bending  is  dependent  upon  the  stiffness  of  the  testing 
frame,  bearing  blocks,  load  cell’,  and  loading  platens,  as  well  as 
the  stiffness  of  the  concrete  test  cube.  Comparison  of  test  results 
in  the  compression-tension  region  with  other  investigators  (Section  VI) 
indicates  that  the  induced  bending  error  may  not  be  ’arge. 

Bending  stresses  were  no  doubt  also  induced  in  Lhe  biaxial  tension 
test  specimens;  however,  their  magnitude  would  be  considerably  less 
than  in  the  compressive-tension  tests.  Bending  stresses  could  be  elim¬ 
inated  by  using  a  testing  machine  designed  such  that  the  loading  axes 
could  be  adjusted  during  a  test. 

4.6  TRIAXIAL  TESTS 

4.6.1  Compression  Tests 

All  of  the  triaxial  tests  were  conducted  using  friction  reducing 

pads. 
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Many  of  the  problems  encountered  during  the  uniaxial  testing  and 
the  biaxial  testing  were  not  applicable  to  triaxial  testing.  All  sides 
were  confined,  therefore,  no  problem  arose  in  preventing  confinement. 

The  effect  of  loading  on  the  different  cube  surfaces  was  not  dis¬ 
covered  until  after  the  triaxial  test  series  had  bt.  2n  completed.  The 
test  records  indicate  that  the  minor  load  was  not  always  applied  to  the 
shaped  surface;  hence,  it  was  concluded  that  some  of  the  scatter  in  the 
test  results  was  due  to  cube  orientation  in  the  testing  machine.  The 
minor  load  should  have  always  been  applied  to  the  shaped  surface. 

The  compressive  capacity  of  the  testing  machine  was  reached  before 
the  test  cube  failed  whenever  the  ratio  of  the  minor  load  (lowest 
absolute  load)  to  the  major  load  (largest  absolute  load)  was  0.3  or 
larger.  At  the  higher  stress  levels,  it  became  difficult  to  operate 
the  pumps  supplying  pressure  to  the  jacks. 

As  in  the  biaxial  tests,  once  the  maximum  stress  had  been  obtained 
on  the  major  axis  (axis  on  which  the  largest  absolute  load  was  applied), 
the  original  loading  proportions  could  not  be  maintained.  It  was  im¬ 
possible  to  decrease  the  load  on  the  minor  axis  once  the  peak  stress 
had  been  reached  on  the  major  axis.  The  test  cube  expanded  rapidly  in 
the  direction  of  the  minor  axis. 

4.6.2  Triaxial  Tension  Tests  and  Combined  Compression-Tension  Tests 

The  elimination  of  bending  stresses  induced  during  the  test  was  the 
greatest  problem  associated  with  triaxial  testing  where  at  least  one  of 
the  loads  was  applied  in  tension.  This  problem  was  not  resolved  in  the 
case  of  triaxial  testing  either. 
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The  triaxial  tension  test  specimens  were  difficult  to  obtain.  The 
tension  heads  had  to  be  glued  to  the  faces  of  the  cubes.  Four  of  the 
tension  heads  could  be  glued  in  place  in  one  operation;  however,  the 
remaining  two  had  to  be  put  in  place  at  a  later  time.  Proper  alignment 
of  the  last  two  tension  heads  was  difficult.  Alignment  was  important 
to  reduce  bending  stresses. 
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SECTION  V 


RESULTS  OF  TESTS 


5.1  INTRODUCTION 

Most  of  the  tests  conducted  in  this  testing  program  were  biaxial 
and  triaxial  compression  tests.  Some  tension-compression  and  tension- 
tension  type  biaxial  tests  were  also  conducted.  In  the  triaxial  series 
of  tests,  some  tension-compression-compression,  tension-tension-compres¬ 
sion  and  tension-tension-tension  type  testing  was  conducted. 

The  strain  data  is  presented  primarily  in  Section  VIII  in  con¬ 
nection  with  the  model  predicted  results  and  in  Appendix  II.  Some  strain 
data  is  also  presented  in  Section  VI  where  comparisons  are  ma>  i  with 
other  published  strain  data. 

5.2  BIAXIAL  STRENGTH  OF  CONCRETE 

All  data  on  the  strength  of  concrete  subjected  to  a  biaxial  state 
of  stress  are  presented  in  normalized  form.  The  results  show  the  ratio 
of  the  maximum  stress  observed  to  the  uniaxial  unconfined  compression 
strength,  o  t  for  various  principal  stress  ratios,  c^/Oj.  ^  value  of  a r 
was  obtained  for  each  batch  of  18  cubes  by  taking  the  average  uniaxial 
unconfined  compression  strength  of  3  cubes  which  were  randomly  selected 
from  the  batch.  Thus,  in  each  batch  of  test  cubes,  15  cubes  were  avail¬ 
able  for  the  various  multiaxial  tests. 

The  relationship  between  principal  stresses,  0^/0^  and  ,  at 

failure  is  shown  in  Figure  13.  In  this  figure,  the  average  biaxial 
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strength  curve  is  shown  based  on  the  biaxial  tests  which  were  conducted. 
Based  on  the  tests  in  the  biaxial  compression  region,  a  strength  of  1.15 
cr  was  observed  for  the  case  of  equal  compression  in  both  principal 
directions.  The  highest  relative  strength  was  observed  at  a  stress  ratio 
of  “2.5  where  c^/o^  “  1.33.  The  curve  is  relatively  flat  in  this 

region  with  a  sharp  rise  from  the  value  of  for  a  very  small  amount  of 
stress  in  the  o ^  direction. 

The  data  points  which  are  the  basis  for  the  biaxial  compression 
curve  shown  in  Figure  13  are  plotted  in  Figure  14  and  tabulated  In 
Table  8.  In  Figure  14  it  can  be  seen  that  most  of  the  data  are  close 
to  the  average  curve.  The  experimental  scatter  which  is  observed  occur¬ 
red  as  a  result  of  the  batch  to  batch  variation  of  the  concrete  and  to 
some  extent  on  the  dimensional  variation  from  cube  to  cube. 

The  biaxial  tension-compression  strength  curve  is  also  shown  in 
Figure  13.  The  data  points  and  curve  to  a  larger  scale  are  plotted  in 
Figure  15.  The  biaxial  tension-tension  strength  curve  is  also  shown  in 
Figures  13  and  15.  The  ratio  of  the  uniaxial  tension  to  the  unconfined 
compression  strength,  Of,  was  approximately  0.11  for  the  batch  of  cubes 
used  in  this  series  of  tests.  In  the  biaxial  tension  region  the 
strength  curve  was  found  to  be  relatively  flat  for  a  range  in  principal 
stress  ratio  =  0  to  1/2.  For  equal  tension  in  both  principal 

directions,  a  strength  of  0.083  a was  observed. 

5.3  TRIAXIAL  STRENGTH  OF  CONCRETE 

The  triaxial  compression  testing  made  up  the  largest  part  of  the 
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triaxial  test  program.  Most  of  the  triaxial  tests  were  conducted  with 
a  ratio  of  the  smallest  compressive  stress,  o^,  equal  to  10%  and  20% 
pf.tjie  largest  compressive  stress,  a^.  For  these  ratios  of  smallest  to 
largest  stress,  varl°us  ratios  of  intermediate  stress,  an  to 

largest  stress,  o^,  were  used.  The  results  of  the  triaxial  compression 
tests  are  presented  in  tabular  form  and  graphical  form. 

In  Table  3  the  average  maximum  values  of  stress  observed  for  tests 
with  various  nominal  stress  ratios  o^/c^/a^  are  shown.  The  normalized 
results  based  on  the  uniaxial  unconfined  compression  strength,  or»  are 
also  tabulated.  The  triaxial  compression  curves  for  *=  0.2  and 
■  O.la-^  are  shown  in  Figure  16.  The  triaxial  compression  strength 
data  are  also  plotted  in  Figure  16.  The  curves  are  a  representation  of 
the  three-dimensional  strength  of  concrete  observed  in  this  test  series. 
The  largest  compressive  stress,  is  plotted  against  the  intermediate 
compressive  stress,  a2*  The  increase  in  the  maximum  principal  stress, 
G^,  is  primarily  a  function  of  the  ratio  of  minimum  principal  stress, 

G^.  to  maximum  principal  stress,  Thus  all  points  with  a  constant 

ratio  o^/a^  will  fall  on  a  curve.  Each  curve  determined  in  this  manner 
will  be  a  line  on  the  three-dimensional  concrete  strength  envelope. 

The  data  plotted  in  Figure  16  are  tabulated  in  Table  9.  The  tri¬ 
axial  compression  strength  data  are  also  listed  in  normalized  form  in 
Table  9.  The  normalized  strength  data  listed  in  Table  9  are  plotted 
in  Figure  17.  The  biaxial  strength  curve,  B  0,  is  also  shown  in 
Figure  17  and  is  seen  as  a  special  case  of  the  triaxial  strength  curves. 
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TABLE  3 


AVERAGE  TRIAXIAL  COMPRESSION  STRENGTH  DATA 


k  -» 

Nominal  Ratio 
°l-02_03 

. V 

(psi) 

- 

(psi) 

'  °3  ‘ 

(psi) 

°l/or 

o0/a 

2  r 

Var 

10-1-1 

10,360 

1,307 

1,120 

3.103 

0.391 

0.335 

10-2-1 

12,650 

2,590 

1,140 

3.569 

0.730 

0.394 

10-3-1 

13,070 

4,070 

1,330 

3.920 

1.220 

0.400 

10-4-1 

10,125 

4,150 

1,000 

2.965 

1.213 

0.293 

10-8-1 

12,800 

10,280 

1,280 

3.647 

2.930 

0.366 

10-10-1 

10,730 

10,660 

1,050 

2.883 

2.871 

0.281 

10-2-2 

19,220 

4,170 

4,000 

5.842 

1.265 

1.215 

10-3-2 

22,130 

6,950 

4,500 

6.641 

2.075 

1.350 

10-5-2 

23,720 

11,820 

4,690 

7.304 

3.640 

1.444 

10-8-2 

24,200 

19,550 

4,810 

7.442 

6.012 

1.479 

10-10-2 

20,660 

20,380 

4,175 

5.703 

5.623 

1.155 

10-10-1/2 

8,283 

8,250 

447 

2.104 

2.097 

0.112 

10-10-3/2 

15,560 

15,500 

2,285 

4.059 

4.045 

0.595 
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Figure  17,  Normalized  Triaxial  Compression  Data 


For  the  case  of  equal  principal  stresses  in  two  directions,  the 
effect  of  a  variation  of  the  third  principal  stress  is  summarized  in 
Table  4.  The  effect  of  a  very  small  amount  of  triaxial  confinement  is 
clearly  seen  here.  For  the  case  of  0^  =  0.05  0  ,  the  triaxial 
strength  is  almost  double  the  biaxial  strength  in  the  direction  of  0^ . 
For  the  case  of  =  0.20  0^,  the  triaxial  strength  is  approximately 
five  times  as  great  as  the  biaxial  strength. 

The  results  of  the  testing  indicate  that  for  a  constant  stress 
ratio,  0  /a^,  the  highest  stress  occurs  approximately  at  a  stress 
ratio  of  =0.5  0^.  In  order  to  establish  this  ratio  more  exactly 
it  would  be  necessary  to  do  more  testing  in  this  region  than  was  done 
on  this  project.  It  would  also  be  desirable  to  use  all  specimens  from 
the  same  batch  of  concrete.  Based  on  the  test  results  obtained,  it 
appears  that  the  stress  increase  in  the  region  =  0.5  0^  is  occurring 
at  a  somewhat  faster  rate  than  the  stress  increase  in  the  region  0^  = 
as  the  third  principal  stress  0^  is  increased. 

Triaxial  tests  in  the  tension-compression  region  were  also  con¬ 
ducted.  The  results  of  the  tests  are  tabulated  in  Table  10  (Appendix 
II).  Since  the  preparation  of  test  specimens  of  this  type  was  very 
difficult,  there  was  no  attempt  made  to  duplicate  results  for  a  given 
principal  stress  ratio.  Individual  tests  were  conducted  for  various 
principal  stress  ratios.  The  stress  ratios  selected  were  varied  to 
cover  as  much  of  this  region  as  possible  and  to  have  at  least  some 
information  over  a  fairly  large  range  of  stress  ratios. 
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TABLE  4 


COMPARISON-  OF-BIAX-IAL- V»TH  TRIAXIAL  STRENQT.il . 


Nominal  Ratio 

°r°2‘a3 

°l/0r 

°2/0r 

a./o 

3  r 

10-10-0 

1.15 

1.14 

0 

10-10-1/2 

2.10 

2.10 

0.112 

10-10-1 

2.88 

2.87 

0.281 

10-10-3/2 

4.06 

4.05 

0.595 

10-10-2 

5.70 

5.62 

1.155 
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The  data  tabulated  in  Table  10  include  the  biaxial  tension-com¬ 
pression  data  alsa. _ These  dqta  are  a  special  case  of  the  triaxial 

tens ion-compress ion  region.  For  biaxial  tension-compression  data,  the 
intermediate  stress,  c^,  is  zero  and  the  principal  stress,  o^,  is  com¬ 
pressive  with  the  principal  stress,  0^,  being  the  tensile  stress. 

For  the  triaxial  tension-compression  data,  the  intermediate  stress, 
o^,  is  either  compressive  or  tensile.  The  data  are  presented  in  norm¬ 
alized  form  with  the  ratios  of  the  principal  stress  to  the  absolute  value 
of  ar  listed.  Thus  each  data  point  has  the  same  sign  as  the  corres¬ 
ponding  principal  stress,  in  each  case. 

The  data  in  Table  10  are  presented  in  groups  which  have  approximately 
the  same  principal  stress  ratio,  a^/o^.  Within  these  groups  it  is 
possible  to  observe  strength  trends  as  the  intermediate  principal  stress, 
02 1  is  varied  from  the  case  O2  11  0^  to  the  other  limiting  case  02  =  0. 

As  in  the  triaxial  compression  region,  the  most  significant  strength 
parameter  is  the  ratio  of  minimum  principal  stress,  0^,  to  maximum 
principal  stress,  o^.  The  variation  of  the  intermediate  principal  stress, 
Oj,  for  a  constant  stress  ratio,  0^/0^,  has  a  secondary  effect. 

As  might  be  expected,  an  increase  in  stress  ratio,  o^/o^,  causes  a 
decrease  in  maximum  compressive  strength  o^.  In  other  words,  a  prin¬ 
cipal  tension  stress  increase  in  one  axis  causes  a  principal  compressive 
stress  decrease  on  another  axis.  For  the  case  of  an  approximately 
equal  stress  ratio,  o^/o^,  an  intermediate  principal  stress  variation 
from  02  *  0  to  ■  0^  causes  a  reduction  in  the  maximum  principal 
stress,  o^. 


A  graphical  representation  of  the  data  tabulated  in  Table  10  is 
shown  in  Figure  18.  Figure  18  is  a  direct  plot  of  the  values  a^/or 
against  the  values  cr_/o  .  The  number  shown  next  to  each  plotted  data 
point  is  the  ratio  c^/o^.  The  biaxial  compression  strength  curve  is 
shown  in  order  to  establish  a  reference  base  since  it  represents  one 
boundary  for  the  region  under  consideration.  The  curves  shown  should 
give  interpolated  values  of  concrete  strength  in  this  region  which  are 
associated  with  the  data  points  through  which  the  curves  pass.  Each 
curve,  therefore,  contains  all  points  on  the  strength  envelope  for  this 
region  which  have  a  constant  ratio  of  a^al‘ 

Table  5  shows  the  triaxial  tension  data  for  the  case  of  three 
approximately  equal  principal  tensile  stresses.  Only  one  test  was 
conducted  in  this  region  since  the  region  is  very  small  compared  to 
other  regions  considered.  It  was  not  expected  that  a  significant  vari¬ 
ation  in  strength  would  occur  for  various  principal  tensile  stress 
ratios.  The  result  in  Table  5  indicates  a  slight  strength  reduction 
compared  to  the  uniaxial  tension  strength  of  approximately  0.11  Of. 

TABLE  5 


TRIAXIAL  TENSION  DATA 


A  final  presentation  of  the  variation  in  triaxial  strength  of 
concrete  is  shown  in  Figure  19.  Two  curves  showing  the  relative  con¬ 
crete  strength  for  various  triaxial  states  of  stress  are  presented. 

The  curves  stiow  how  the  ratio  of  minimum  compressive  stress  to 
maximum  compressive  stress  affect  the  maximum  strength  obtained 
during  a  test.  The  two  cases  included  are  the  cases  where  the  inter¬ 
mediate  compressive  stress  and  a ^  =  0.5  a^.  For  both  cases, 

the  curves  show  the  decrease  in  for  a  decrease  in  the  ratio  o^/o^ 
in  the  triaxial  compression  region.  In  the  triaxial  tension-compres¬ 
sion  region,  there  is  a  continued  decrease  in  as  the  ratio  of  ten¬ 
sile  stress  to  compressive  stress  increases.  The  points  on  the 
curve  where  =  0  are  points  on  the  biaxial  compression  curve. 

5.4  FAILURE  MODES 

The  crack  patterns  observed  in  the  specimens  after  failure  (after 
the  ultimate  strength  had  been  obtained)  were  similar  to  those  reported 
by  previous  investigators. 

Numerous  cracks  in  a  direction  parallel  to  the  applied  loads  were 
observed  in  cubes  subjected  to  uniaxial  compression.  Several  examples 
of  failures  of  this  type  are  shown  in  Figure  20.  The  top  side  of  the 
cube,  as  seen  in  Figure  20,  was  the  loaded  surface.  The  crack  pattern 
on  the  loaded  surface  is  random. 

In  cubes  subjected  to  equal  biaxial  compression,  the  cracking 
occurred  on  planes  parallel  to  the  unloaded  surface.  This  type  of 
crack  pattern  develops  since  the  cube  is  expanding  uniformly  in  the 
direction  in  which  there  is  no  confinement.  Failure  occurs  when  this 
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Figure  19.  Triaxial  Concrete  Strength 
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Figure  20.  Uniaxial  Type  Failures 
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expansion  becomes  excessive.  In  Figure  21  cubes  with  this  type  of 
failure  pattern  are  shown.  The  top  surface  of  the  cubes  was  the 
free  surface.  The  cracks  on  the  loaded  surfaces  are  numerous  and 
in  general  parallel  to  each  other.  Cubes  loaded  in  unequal  biaxial 
compression  are  shown  in  Figure  22.  The  cracking  observed  here  was 
somewhat  different  from  the  cracking  observed  in  the  cubes  subjected 
to  equal  biaxial  compression.  As  the  stress  state  changes  from 
close  to  equal  to  unequal  biaxial  compression,  the  cracking  changes 
from  numerous  parallel  failure  planes  to  a  few  major  nonparallel 
failure  planes.  The  angle  between  the  failure  planes  increases  as 
the  stresses  become  more  unequal.  The  cracks  on  the  surface  with  the 
higher  stress  are  essentially  parallel,  with  the  nonparallel  cracks 
occurring  on  the  surface  with  the  lower  stress. 

Cubes  which  were  subjected  to  high  triaxial  compression  are 
shown  in  Figure  23.  The  cubes  were  unloaded  carefully  after  the  max¬ 
imum  strength  had  been  attained  for  the  given  loading  ratios.  Large 
dimensional  changes  have  occurred.  Numerous  random  microcracks  exist 
in  all  directions.  The  cubes  will  crumble  if  they  are  not  handled 
carefully.  The  observed  failure  modes  for  cubes  subjected  to  small 
triaxial  compressive  stresses  fall  into  three  categories.  For  the 
case  where  the  two  smaller  principal  stresses  are  equal,  the  failure 
is  as  shown  in  Figure  20.  For  the  case  where  the  two  larger  principal 
stresses  are  equal,  the  failure  is  of  the  type  shown  in  Figure  21. 

For  the  case  where  the  three  principal  stresses  are  all  different, 
the  failure  was  by  splitting  in  the  direction  of  the  minor  principal 
stress  similar  to  that  shown  in  Figure  22. 
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Figure  21.  Equal  Biaxial  Type  Failures 
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Figure  22.  Unequal  Biaxial  Type  Failures 


triaxial 


Figure  23.  Triaxial  Type  Failures 
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Failure  in  the  cubes  subjected  to  uniaxial  tension  was  due  to 
a  major  crack  through  the  specimen  which  was  essentially  parallel  to 
the  loaded  surfaces.  The  orientation  of  the  failure  plane  in  biaxial 
tension  varied  from  0-45  degrees  with  the  horizontal  as  the  tensile 
stresses  varied  from  the  uniaxial  case  to  the  equal  biaxial  stress 
case.  In  the  "unos  subjected  to  triaxial  tension,  the  failure  was 
due  to  a  single  major  crack  with  no  preferred  direction. 

5.5  CUBE  VERSUS  CYLINDER  STRENGTHS 

A  tabulation  of  the  unconfined  compression  strength,  a  ,  and  the 
7-  and  28-day  cylinder  strengths  is  given  in  Table  6.  The  cylinder 
tests  were  used  as  initial  control  on  each  batch.  Once  a  strength 
curve  was  determined,  it  was  possible  to  predict  whether  the  batch 
would  conform  to  the  specified  strength  based  on  the  7-day  tests. 

The  28-day  strengths  were  taken  4  weeks  after  the  casting  date  while 
the  values  of  were  taken  at  the  test  date.  It  was  initially 
assumed  that  there  would  not  be  an  appreciable  strength  increase  in 
the  concrete  after  28  days  since  all  specimens  were  removed  from  the 
curing  room  and  stored  in  a  dry  condition  until  tested.  Tests  on 
two  control  cylinders  which  were  broken  at  the  time  that  the  uniaxial 
cube  strengths  were  determined  indicated  a  14%  to  20%  increase  in  the 
strength  of  the  batch.  Since  the  cubes  were  tested  approximately 
four  months  after  the  cylinders  in  each  case,  a  direct  strength  com¬ 
parison  cannot  be  obtained.  Another  factor  which  makes  a  direct 
strength  comparison  difficult  is  that  the  cube  strength  is  based  ori 
a  cube  with  beveled  edges  In  each  case.  This  gives  an  apparent  cube 
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TABLE  6 


UNIAXIAL  TEST  RESULTS 


6"  x  12"  Cylinder 


Cube 

Strength 

Date 

Tested 

Date 

Cast 

0 

r 

(psi) 

7 -day 
(psi) 

28-day 

(psi) 

5-20-71 

2-3-71 

3,620 

2,953 

4,314 

5-18-71 

2-2-71 

3,412 

3,200 

4,173 

6-29-71 

2-9-71 

3,756 

3,236 

4,420 

6-7-71 

2-8-71 

3,930 

3,483 

4,527 

6-16-71 

2-4-71 

3,340 

3,324 

4,314 

7-2-71 

2-11-71 

3,333 

2,758 

3,996 

7-2-71 

3-2-71 

3,327 

2,599 

4,029 

7-12-71 

2-16-71 

3,500 

3,183 

4,332 

7-20-71 

3-8-71 

3,307 

3,271 

4,827 

7-28-71 

3-11-71 

4,337 

3,094 

4,810 

8-3-71 

3-10-71 

3,161 

3,359 

4,810 

8-10-71 

4-7-71 

2,597 

3,554 

4,845 

8-19-71 

4-12-71 

2,977 

3,324 

4,951 
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strength  which  is  too  low.  It  is  felt  that  a  10%  to  12%  increase 
in  the  cube  strength  would  give  a  better  comparison  basis. 

In  general  the  cube  strength  falls  between  the  7-day  and  28-day 
cylinder  strengths.  For  the  last  three  batches  shown  in  Table  6, 
the  cube  strengths  were  lower  than  the  7-day  cylinder  strengths. 

This  inconsistency  could  be  due  to  any  number  of  factors  such  as 
difference  in  compaction,  segregation  of  material  or  curing  method. 
Since  the  cube  strength  in  each  case  was  used  as  the  strength  basis, 
the  inconsistency  noted  here  is  not  of  great  importance. 

5.6  PLATEN  SHEARING  EFFECT 

It  was  necessary  to  use  a  loading  platen  which  was  smaller  than 
the  cube  in  order  to  allow  for  deformation  of  the  specimen  during  a 
test.  In  order  to  obtain  a  better  stress  distribution,  the  cubes  were 
beveled  along  the  edges.  For  the  4000  psi  concrete  mix  used  in  this 
testing  program  and  for  cubes  tested  using  friction  reduction  pads, 
a  shearing  effect  along  the  edge  of  the  platen  was  not  observed  to  be 
a  significant  problem  in  this  test  program.  If  shearing  was  signif¬ 
icant,  thin  plates  of  concrete  would  have  broken  off  on  all  sides 
during  a  uniaxial  test.  Typical  uniaxial  type  failures  are  shown  in 
Figure  20.  The  type  of  failure  described  above  did  not  occur. 
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SECTION  VI 


COMPARISON  OF  DATA  WITH  OTHER  INVESTIGATORS 

6.1  INTRODUCTION 

A  comparison  of  published  data  with  the  data  obtained  from  this 
testing  program  can  be  accomplished  generally  for  strength  comparisons 
only.  In  only  one  case,  Hilsdorf  (4),  are  there  strains  published  that 
can  be  compared.  This  is  for  the  biaxial  case  only.  Mills  (25)  has  a 
few  stress-strain  curves  for  the  triaxial  case,  but  not  enough  to  be 
able  to  compare  completely. 

In  comparing  the  ultimate  strengths  from  the  various  multiloading 
test  results  published,  difficulty  is  encountered  in  that  different 
authors  present  their  results  differently.  Some  use  normalized  results 
whereas  others  do  not.  Still  another  complication  arises  since  different 
types  of  test  specimens  were  used  and  the  stresses  or  loads  were  applied 
differently.  In  some  cases  cylindrical  specimens  were  used,  in  other 
cases  slabs  were  used,  and  in  others  cubes  were  used.  Hollow  as  well 
as  solid  cylinders  were  used.  Hollow  cylinders  were  generally  loaded 
axially  along  with  a  torsion  or  fluid  pressure.  Solid  cylinuers  were 
generally  tested  in  a  solid  type  triaxial  testing  machine.  Slabs  were 
tested  biaxially  only.  Only  those  using  cube  specimens  could  indepen- 
ently  control  the  applied  load  in  three  directions.  Usually  only  the 
ultimate  strengths  were  reported.  Measuring  of  strains  was  difficult 
or  impossible  in  many  of  the  different  types  of  test  setups. 
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6.2  BIAXIAL  COMPRESSION  STRENGTH  COMPARISON 


A  comparison  of  the  biaxial  strength  of  concrete  determined  from 
this  test  program  with  the  results  obtained  by  other  investigators  has 
been  made.  In  Figure  24,  the  results  of  the  biaxial  compression  investi¬ 
gators  were  determined  as  accurately  as  possible  from  the  available  re¬ 
ports.  In  some  cases  the  published  data  was  very  limited  but  the 
curves  in  general  show  the  reported  trends. 

Rosenthal-Glucklich  (22),  Campbell-Alien  (14)  and  Bellamy  (11) 
used  hollow  cylinders  subjected  to  axial  compression  and  external  hydrau¬ 
lic  pressure.  For  the  case  of  equal  biaxial  compression,  the  range  of 
values  for  cr^/c^  was  2.24  to  2.69.  A  triaxial  state  of  stress  is  sus¬ 
pected  in  this  type  of  test  and  the  determined  strength  increase  is 
considered  to  be  too  high. 

Iyengar  et  al.  (5)  used  4-in.  and  6-in.  cubes  with  solid  bearing 
platens  and  no  friction  reduction  attempted  between  the  specimen  and 
the  platens.  Using  this  type  of  testing  procedure,  the  friction  between 
specimen  and  platen  has  a  confining  effect  on  the  concrete  and  will  re¬ 
sult  in  an  overestimate  of  strength  increase.  For  the  case  of  equal 
biaxial  compression,  these  tests  resulted  in  a  ratio  c^/o^.  *  3.6.  This 
estimate  of  the  biaxial  compression  strength  is  not  as  good  as  that 
obtained  using  hollow  cylinder  test  specimens. 

Tests  using  square  slab-type  specimens  subjected  to  in-plane  load¬ 
ing  were  conducted  by  Vile  (26)  and  Robinson  (3).  Both  investigators 
used  10  x  10  x  4-inch  specimens  and  solid  bearing  platens  with  a  special 
concrete  curing  compound  used  as  packing  to  reduce  the  friction  effects 
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Figure  24.  Biaxial  Compressive  Strength  Comparison 


between  the  specimen  and  the  platens.  For  the  case  of  equal  biaxial 
compression,  a  ratio  o^/o^  of  1.10  was  obtained  by  Robinson,  and  1.25 
to  1.45  was  obtained  by  Vile.  Some  of  the  difference  in  results  ob¬ 
served  by  these  two  investigators  was  probably  due  to  the  different 
concrete  mixes  which  were  used. 

Weigler  and  Becker  (27)  used  10  x  10  x  2.5-cm  slabs  with  solid 
bearing  platens  and  no  friction  reduction  material.  Their  investiga¬ 
tions  yielded  a  ratio  of  a^/or  of  1.10  to  1.20  for  the  case  of  equal 
biaxial  compression. 

Another  investigation  using  a  slab  type  test  specimen  was  conducted 
by  Kupfer,  Hilsdorf  and  Rusch  (4).  A  20  x  20  x  5-cm  specimen  was  used. 
This  research  made  use  of  a  brush  type  bearing  platen  in  an  attempt  to 
load  the  test  specimen  without  end  restraint.  A  ratio  of  o^/o^  *  1.16 
was  found  for  the  case  of  equal  biaxial  compression. 

Mills  and  Zimmerman  (21)  used  a  2  1/4-inch  cube  for  a  test  speci¬ 
men.  In  order  to  eliminate  end  constraint,  a  system  using  two  polyethy¬ 
lene  pads  with  grease  between  the  pads  was  used.  The  biaxial  tests,  in 
this  investigation,  vere  conducted  using  a  sequential  loading.  The 
case  of  equal  biaxial  loading  was  not  quite  achieved  using  this  loading 
method.  Three  different  mixes  were  used  with  average  values  of  of 
3,340,  3,910  and  5,235  psi,  respectively.  The  biaxial  tests  for  the 
two  higher  strength  mires  were  too  few  to  achieve  reliable  results.  For 
the  lower  strength  mix,  a  value  of  O^/o^  for  the  case  of  equal  biaxial 
compression  is  approximately  1.24. 
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From  the  previous  discussion  it  can  be  concluded  that  all  of  the 
previous  biaxial  compression  studies  did  not  use  a  test  set-up  which 
developed  unconfined  biaxial  state  of  stress.  The  hollow  cylinder  and 
cube  specimens  with  solid  bearing  platens  and  no  friction  reduction 
between  the  concrete  cube  and  the  bearing  platen  fall  into  this  class. 

The  investigo cions  using  slab  type  specimens  with  solid  and  brush 
type  bearing  platens  seem  to  develop  a  reasonably  unconfined  biaxial 
stress  state.  The  results  of  thc?e  investigations  are  comparable  to 
the  results  obtained  from  this  testing  program  using  a  cube  type  test 
specimen  with  greased  pads  to  eliminate  platen  restraint.  The  biaxial 
compression  strength  curve'  obtained  using  the  slab  type  specimens  are 
all  within  10%  of  the  results  obtained  during  the  present  investigation. 
This  appears  to  be  a  reasonable  variation  since  the  concrete  mixes 
used  in  all  the  investigations  were  different.  In  fact,  during  the 
present  investigation,  a  different  batch  of  concrete  using  the  same 
mix  was  found  to  vary  almost  as  much. 

The  results  obtained  using  the  brush  type  bearing  platen  with 
slab  type  specimen  compare  best  to  the  results  obtained  in  this  test 
program.  This  is  not  unexpected  since  of  all  other  systems  studied 
it  probably  achieves  the  most  unconfined  biaxial  compression  state  of 
stress. 

6.3  BIAXIAL  TENSION-COMPRESSION  STRENGTH  COMPARISON 

Previous  investigations  of  the  biaxial  tension-compression  strength 
of  concrete  were  conducted  on  two  different  types  of  specimens. 
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Hollow  cylinders  loaded  axially  along  with  a  fluid  pressure  was 
used  by  Rosenthal-Glucklich  (22),  McHenry-Karni  (8)  and  Tsuboi-Suenaga 
(18). 

Kupfer,  Hilsdorf,  and  Rusch  (4)  used  a  slab  type  specimen  20  x 
20  x  5-cm  subjected  to  an  in-plane  loading.  Brush  type  loading  platens 
were  used  to  reduce  end  restraint  with  the  specimen  glued  to  the  load¬ 
ing  platen  for  application  of  the  tensile  stress. 

A  typical  strength  curve  from  each  of  the  above  test  programs 
is  compared  with  the  results  obtained  from  the  present  investigation 
in  Figure  25.  The  results  in  Figure  25  are  in  normalized  form  in 
order  to  make  a  direct  comparison.  The  stresses  are  normalized  by 
taking  the  ratio  of  each  stress  to  the  unconflned  uniaxial  compressive 
stress  of  concrete,  o^.  In  all  of  the  above  test  programs,  tests  on 
different  concrete  mixes  were  conducted.  The  curves  using  different 
mixes  with  the  same  test  set-up  have  the  same  general  shape.  Another 
trend  which  was  observed  is  that  the  ratio  of  uniaxial  tension  to  uni¬ 
axial  compression  decreases  with  an  increase  in  compression  strength 
for  various  mixes.  The  range  of  values  observed  during  the  previous 
test  programs  for  the  ratio  of  uniaxial  tension  to  the  absolute  value 
of  uniaxial  compression  was  0.07  to  0.13.  This  compares  to  a  value 
of  0.11  obtained  during  this  investigation.  This  value  appears  to  be 
reasonable  since  the  compressive  strength  of  the  mix  used  was  within 
the  range  of  compressive  strengths  used  by  the  previous  investigators. 
In  the  region  where  the  ratio  of  compression  to  tension  is  high,  the 
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Figure  25.  Biaxial  Tension  and  Tension-Compression  Strength  Comparison 


strength  curve  obtained  is  lower  than  other  reported  results.  However, 
in  general  there  is  not  too  much  difference  between  these  results  and 
those  of  Rosenthal  (22) .  This  is  encouraging  since  of  all  the  cylinder 
type  testing  the  diameter  ratio  in  that  test  program  was  closest  to 
thin  wall  assumptions  which  were  used  in  all  the  hollow  cylinder  test 
programs. 

6.4  BIAXIAL  TENSION  STRENGTH  COMPARISONS 

Biaxial  tension  strength  testing  was  previously  conducted  by  Kupfer, 
Hilsdorf  and  Rusch  (4) .  They  used  the  same  test  set-up  here  as  for  the 
tension-compression  testing  except  that  the  specimen  was  glued  to  the 
platens  on  both  sides  in  order  to  apply  a  tension  stress  in  orthogonal 
directions.  A  very  flat  biaxial  tension  curve  was  reported  using  this 
testing  method.  Figure  25  shows  the  biaxial  tension  comparison  curves. 

It  was  reported  that  the  biaxial  tensile  strength  is  approximately 
equal  to  Its  uniaxial  strength. 

Rosenthal  and  Glucklich  (22)  conducted  biaxial  tension  tests  by 
loading  hollow  cylinders  with  axial  tension  and  internal  fluid  pressure. 
The  resulting  biaxial  tension  strength  shows  a  small  decrease  in  the 
biaxial  tension  strength  compared  to  its  uniaxial  tension  strength. 

6.5  TRIAXIAL  STRENGTH  COMPARISONS 

A  comparison  of  the  t.iaxial  strength  of  concrete  can  be  made  in 
the  triaxial  compression  region  only.  There  have  been  no  other  pub¬ 
lished  test  results  for  the  triaxial  tension-compression  region.  In 
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the  triaxial  compression  region,  it  is  difficult  to  get  a  direct 
strength  comparison  since  there  are  so  many  different  loading  combina¬ 
tions  which  are  possible.  Most  of  the  testing  which  was  performed  by 
others  used  loading  combinations  which  were  different  from  those  used 
in  this  test  program.  There  was  also  a  difference  in  types  of  speci¬ 
mens  used,  concrete  mixtures,  end  conditions  and  the  uniaxial  uncon¬ 
fined  compression  strength.  Any  of  these  factors  couJd  cause  some 
difference  in  the  results  obtained  by  the  other  investigators. 

The  specimens  used  by  other  investigators  include  cubes,  solid 
cylinders  and  hollow  cylinders.  A  summary  of  other  investigations 
including  specimen  type,  maximum  aggregate  size  and  end  conditions  is 
shown  in  Table  7. 


TABLE  7 


SUMMARY  OF  PREVIOUS  TRIAXIAL  COMPRESSION  INVESTIGATORS 


Investigators 


Specimen  Type 


Maximum  Size 
Aggregate 


End 

Condition 


Balmer  (13) 

Bellamy  (11) 

Gardner  (2) 

Krishna swamy  (1) 
Krishnaswamy  (1) 
Krishnaswamy  (1) 

Mills  &  Zimmerman  (21) 


6"xl2"  Solid  Cylinders 
6"xl2"  Hollow  Cylinders 
3"x6"  Solid  Cylinders 
4"  Cube 
4"  Cube 
4"  Cube 
2  1/4"  Cube 


1  1/2  in. 
Sand 
3/4  in. 
3/4  in. 
3/4  in. 
3/4  in. 
3/8  in. 


Steel  Plates 
Not  Specified 
Not  Specified 
No  Capping 
Plastic  Sheet 
Plastic  &  Grease 
Plastic  &  Grease 


The  data  obtained  from  these  testing  programs  are  limited  to 


a  very  small  portion  of  the  triaxial  compression  region.  The 
results  of  tests  using  cylinders  is  shown  in  Figure  26. 


80 


81 


Figure  26.  Triaxial  Compression  Strength  Comparisons 


The  triaxial  compression  curve  obtained  during  this  investigation  is 
also  shown  in  Figure  26  so  that  a  comparison  can  be  shown. 

For  all  cylinder  tests,  a ^  =  O ^  due  to  nature  of  the  test  pro¬ 
gram.  Since  the  maximum  axial  strength,  a^,  increases  very  rapidly 
for  small  confining  pressures,  0 ^  and  0^,  the  system  becomes  self- 
limiting  to  a  relatively  small  region.  The  highest  ratio  shown  in 
Figure  26  is  0^/0^  *  0.27.  Based  on  the  triaxial  compression  curves 
obtained  during  the  present  investigation,  the  increase  in  the  maximum 
compressive  strength,  0^  reported  by  Bellamy  (11)  appears  to  be  too 
low  for  all  reported  results.  The  results  reported  by  Gardner  are 
somewhat  low  for  the  ratio  o^/o^  ■  0.12  but  seem  to  compare  quite  well 
for  the  other  two  data  points.  The  results  reported  by  Balmer  (13) 
are  a  bit  high  for  all  reported  data. 

The  comparison  of  results  for  tests  which  used  a  cube  test  speci¬ 
men  are  shown  in  Figure  27.  Using  a  cube  test  specimen  all  three  prin¬ 
cipal  stresses  can  be  varied.  It  is  thus  possible  to  obtain  test  data 
over  a  larger  portion  of  the  triaxial  compression  region.  Test  results 
reported  by  other  Investigators  is  limited  to  the  lower  portion  of  the 
region  due  to  loading  limitations  of  the  test  set-up  which  was  used. 

The  results  are  limited  to  a  principal  stress  ratio,  0^/0^*  approxi¬ 
mately  0.13. 

Based  on  the  triaxial  compression  curves  obtained  during  the  pres¬ 
ent  investigation,  a  comparison  can  be  made  with  the  results  obtained 
by  Mills  and  Zimmerman  (21)  for  a  stress  ratio,  0^/0^  of  approximately 
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Figure  27.  Triaxial  Compression  Strength  Comparisons 


0.10.  The  data  reported  by  them  seem  to  fall  on  two  curves.  One 
curve  is  quite  a  bit  lower  than  the  other.  The  upper  curve  shown  in 
Figure  27  compares  quite  well  with  the  results  of  the  present  investi¬ 
gation. 

The  results  reported  by  Krishnaswamy  (1)  were  based  on  tests  with 
end  conditions  which  resulted  in  various  degrees  of  constraint  of  the 
cube  by  the  loading  platens.  The  various  conditions  were  two  plastic 
sheets  and  grease,  one  plastic  sheet  and  no  capping.  His  reported 
results  of  tests  using  a  stress  ratio,  0^/0^,  °f  approximately  0.10 
are  shown  in  Figure  27.  Using  either  one  plastic  sheet  or  two  plastic 
sheets  with  grease  between  the  sheets  did  not  produce  a  significant 
difference  in  results.  The  results  reported  indicate  a  much  lower 
strength  Increase  for  this  type  of  test  than  for  the  tests  where  no 
capping  was  used.  For  the  case  where  no  capping  was  used  the  strength 
increase  Indicated  was  somewhat  less  than  that  found  during  the  present 
investigation. 

A  comparison  of  the  strength  increase  for  the  triaxial  stress  case 
where  0^  =  0^  was  attempted.  There  was  not  a  sufficient  amount  of  pub¬ 
lished  data  to  make  such  a  comparison.  It  was  also  not  possible  to 
make  any  comparisons  of  triaxial  compression  strength  for  the  curve 
where  “  0.20  0^  in  the  region  between  o ^  *  0.20  0^  and  0 ^  ■  0^. 

It  was  also  not  possible  to  make  any  comparisons  in  the  triaxial 
tension  and  triaxial  tension-compression  region  since  no  data  of  this 
type  has  been  reported  by  other  investigators. 
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6.6  COMPARISON  OF  STRESS-STRAIN  CURVES 


A  comparison  of  stress-strain  curves  with  Hilsdcrf  et  al.  (4) 
was  possible.  Different  strength  concretes  were  involved,  hence, 
stresses  are  plotted  on  a  comparable  basis.  These  curves  are  shown 
in  Figure  28.  For  the  uniaxial  compression  case,  the  stress-strain 
curves  from  these  tests  and  those  of  Hilsdorf  et  al.  (4)  were  quite 
similar  in  the  lower  stress  range  but  departed  at  the  higher  stress 
levels.  This  can  be  explained  by  the  time  element  involved  in  the 
testing  procedure.  Hilsdorf  et  al.  required  about  20  minutes  for  a 
test;  whereas,  these  required  about  one  minute.  The  larger  strains 
obtained  by  Hilsdorf  et  al.  can  be  explained  by  creep  effects  which 
are  more  pronounced  at  the  higher  stress  levels.  The  lateral  strains 
also  compare  favorably  at  the  lower  stress  levels,  but  again  depart 
at  the  higher  stress  levels.  In  the  case  of  the  lateral  strains,  the 
results  of  Hilsdorf  are  smaller  than  these.  This  could  possibly  be 
explained  by  the  different  types  of  test  specimens  used.  Hilsdorf  used 
square  slabs  whose  thickness  was  much  smaller  than  the  other  dimensions. 
There  were  probably  fewer  cracks  forming  near  failure,  hence  a  smaller 
apparent  lateraJ  strain  was  detected.  In  the  cube  specimens,  the 
lateral  strain  in  one  unloaded  direction  always  exceeded  the  lateral  strain 
in  the  other  unloaded  direction.  Failure  first  occurred  in  the  direction 
with  the  greatest  lateral  strain... 

In  the  case  of  biaxial  compression  with  the  two  applied  stresses 
being  equal,  again  the  results  of  Hilsdorf  and  the  author's  compare 
favorably  at  the  lower  stress  levels  and  depart  at  the  higher  stress 
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levels.  The  same  trends  are  observed  as  in  the  case  of  uniaxial  tests 


for  both  lateral  and  axial  strains,  however,  the  difference  is  more 
pronounced.  The  lateral  strain  in  these  tests  departs  even  at  the 
lower  loads.  This  is  believed  to  be  due  to  the  different  types  of 
test  specimen  used. 

In  Figure  29  a  comparison  is  made  for  the  case  of  uniaxial  ten¬ 
sion.  Again  the  results  are  quite  comparable.  The  stress-strain 
curve  lias  the  same  shape  with  slightly  smaller  strains  obtained  by 
these.  Figure  29  also  contains  typical  biaxial  tension-compression 
stress-strain  curves.  There  is  not  much  difference  in  the  initial 
portion  of  the  compression  axis  curve.  On  the  tension  axis  however, 
the  curves  depart  at  a  low  stress  level.  A  comparison  with  Hilsdorf 
et  al.  (4)  is  difficult  since  different  ratios  of  tension  to  compres¬ 
sion  were  used. 

For  the  case  of  biaxial  tension,  the  stress-strain  curves  ob¬ 
tained  during  the  present  investigation  are  shown  in  Figure  30.  The 
strain  on  each  loaded  axis  for  the  stress  ratios  cr^/c^  "  *  ant*  ^ 
are  shown. 


6.7  OTHER  COMPARISONS 

Nearly  all  of  the  published  results,  whether  biaxial  or  triaxlal 

could  be  compared  simultaneously  by  plotting  the  mean  normal  stress 

(a  )  versus  the  octahedral  shear  stress  (o  )  at  failure  where: 
o  o 


°o  mJ(0l  + 


a„  + 


°3) 


and 
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Figure  30.  Biaxial  Tension  Stress-Strain  Curves 


T 

O 


1 

y  t (ox  -  a2)2  +  (o2  -  o3)2  +  (o1  -  o3)2]2 


This  was  done  and  some  of  the  published  results  are  shown  in  Figure 
31.  A  rather  large  scatter  is  observed  and  the  results  of  the  pres¬ 
ent  testing  program  fall  within  the  range  of  all  published  results. 
Only  test  results  designated  as  Type  I  (a^  >  *  a^)  and  Type  II 

(o^  m  °2  >  are  comPared«  Tests  with  the  intermediate  principal 
stress,  a2,  different  from  0^  or  will  fall  between  the  two  curves 

shown  in  Figure  31. 

A  comparison  of  mean  normal  stress,  o^,  versus  volumetric  strain, 
Av/v,  is  possible.  Curves  of  this  type  have  been  plotted  in  Figures 
32  and  33.  In  Figure  32  the  results  of  some  typical  biaxial  compres¬ 
sion  tests  are  shown.  In  Figure  33  the  resultB  of  some  triaxial  com¬ 
pression  tests  are  shown  along  with  a  typical  biaxial  and  uniaxial 
test  for  comparison. 
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Figure  32.  Mean  Normal  Stress  versus  Volumetric  Strain  for  Biaxial  Compression 


(ksi) 
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Figure  33.  Mean  Normal  Stress  versus  Volumetric  Strain  for  Triaxial  Compression 


SECTION  VII 


THE  MODEL 


7.1  INTRODUCTION 

Several  attempts  were  made  at  formulating  constitutive  relationships 
for  concrete  subjected  to  a  general  loading  during  this  investigation. 

The  earliest  attempt  utilized  the  linear  constitutive  relationships  of 
the  theory  of  elasticity  with  nonlinearities  introduced  by  considering 
the  modulus  of  elasticity  and  Poisson's  ratio  as  functions  of  the  stresses 
and  strains  instead  of  as  constants.  The  slightest  degree  of  success  was 
never  achieved  using  this  approach;  hence,  it  was  abandoned. 

A  model  to  simulate  concrete  behavior  under  triaxial  loads  was  then 
considered.  The  model  consisted  of  nodal  points  located  at  the  corners 
of  a  cube  along  with  "bars"  or  "springs"  connecting  the  nodal  points.  Any 
number  of  these  cubical  model  elements  could  be  stacked  to  form  a  more 
sophisticated  model.  This  model  could  handle  shear  stresses  and  strains 
in  addition  to  the  normal  stresses  and  strains.  It  could  be  used  to  simu¬ 
late  a  cubical  test  specimen  and  would  be  useful  for  studying  friction  ef¬ 
fects  and  failure  modes  provided  that  strains  were  determined  from  the 
stresses.  Considering  the  simpliest  version  of  this  model,  there  are  eight 
nodes  with  three  possible  displacements  at  each  node;  hence,  a  total  of 
twenty-four  possible  displacements.  A  maximum  of  six  strain  values  exist  at  a 
point.  Conversion  of  six  known  strains  to  twenty-four  displacements  was  con¬ 
sidered  an  impossible  task.  It  was  concluded  that  this  model  could  not  be 
used  to  determine  stresses  whenever  strains  were  given;  hence,  it  was  aban¬ 
doned. 
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The  proposed  model  for  simulating  concrete  behavior  also  consists  of 
nodal  points  connected  by  "bars."  The  nodal  points  are  located  at  the 
points  of  an  octahedron.  The  arrangement  of  the  nodal  points  and  bars  is 
shown  in  Figure  34.  This  arrangement  was  patterned  after  the  model  used 
by  Anson  (23)  in  his  study  of  the  failure  mechanism  for  concrete.  The 
behavior  of  Anson's  model  and  the  one  described  herein  is  quite  different. 
This  model  cannot  handle  shear  stresses  or  strains;  hence,  it  will  be 
developed  with  reference  to  principal  stress  coordinates. 

7 . 2  DEVELOPMENT  OF  THE  MODEL 

The  nomenclature  used  in  conjunction  with  the  nodal  points,  bars,  axes, 
and  dimensions  for  the  model  are  indicated  in  Figure  34.  The  bar  stiff¬ 
nesses  are  shown  below.  The  bars  are  indicated  by  the  nodal  points  at  the 
ends  of  the  bars. 


Bars  (Between  Nodal  Points) 
1-1 
1-2 

1- 3 

2- 2 

2- 3 

3- 3 


Stiffness 


The  model  was  developed  using  forces  and  displacements  instead  of 
stresses  and  strains.  The  forces  and  displacements  are  related  to  the 
stresses  and  strains,  respectively,  as  follows: 
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(7.1) 


where  the  Fs  are  forces,  the  Xs  are  displacements,  u,  v,  and  w  are  the 
model  dimensions,  the  Os  are  stresses,  and  the  £s  are  strains. 

The  mathematical  solution  was  formulated  using  the  stiffness  method 
of  structural  analysis. 

The  relations  between  the  forces  and  displacements  in  the  individual 
members  (bars)  of  the  model  are  expressed  in  matrix  notation  as 


{f}15xl  tkl15xl5  {x}15xl 


(7.3) 


where 

{ f }  is  the  bar  column  matrix 
{x}  is  the  bar  displacement  column  matrix,  and 
[k]  is  the  diagonal  element  stiffness  matrix. 
Matrix  [k]  is  shown  below. 
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The  relations  between  the  bar  displacements  and  the  nodal  displace¬ 
ments  are  expressed  as 


(7.5) 


where 

(x)  is  the  nodal  displacement  column  matrix,  and 
[A]  is  the  rectangular  compatibility  matrix. 
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The  compatibility  matrix  in  terms  of  the  mode]  dimensions  u,  v,  and 
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The  relations  between  the  nodal  forces  and  the  bar  forces  is  expressed 


as 


{F}3xl  =  [ii]3xl5 


{f}15xl  =  lA 


'3x15 


{f} 


15x1 


(7.8) 


where 


{ F }  is  the  nodal  force  column  matrix 

[ B ]  is  the  rectangular  equilibrium  matrix,  and 
T 

[A  ]  is  the  transpose  of  A,  the  capability  matrix. 

The  equilibrium  matrix  is  easily  obtained  by  taking  the  transpose  of 
the  compatibility  matrix. 

Substituting  Equation  7.5  into  Equation  7.3  yields 

{f}  =  [k][A]{X>  (7.9) 

Substituting  Equation  7.9  into  Equation  7.8  yields 

{F}  =  [AT][k][A]{x>  =  [K] (X>  (7.10) 


where 


[K]  is  the  square  3x3  model  stiffness  matrix. 
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Now 


[K]  =  [AT)(k][A] 


(7.11) 


hence 


[K]  = 


4kl02  +  4k6<F)2  +  4k2  4kl  ~2 

a 


4k 


uw 
6b2 


4k 


uv 
Xb2 


4ki  (“) 2  +  4k3£)2  +  4k4 


vw 

4R3  ~~1 

c 


uw 

6 12 


4k 


vw 
3  c2 


4k6(E)2  +  4k3<?>2  +  4k5 


(7. 1.?) 


For  the  model  to  yield  isotropic  results,  the  following  conditions 
were  necessary: 

1.  The  model  dimensions  u,  v,  and  w,  had  to  be  equal  to  each  other. 

2.  The  bar  stiffnesses  kp  k^,  and  k^  had  to  be  equal  to  each  other 
ini tially. 

3.  The  bar  stiffnesses  k^,  k^,  and  k,.  had  to  be  equal  to  each  other 
initially. 

The  stiffnesses  were  selected  such  that  the  stress-strain  results 
yielded  by  the  model  correspond  to  the  stress-strain  results  obtained 
from  tests.  This  is  explained  in  more  detail  later.  The  only  requirements 
on  the  model  dimensions  u,  v,  and  w  were  that  they  be  equal.  This  was 
required  for  isotropic  results  as  mentioned  earlier.  An  arbitrary  and 
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convenient  value  of  10  units  was  therefore  assigned  to  each  of  the  model 
dimensions.  As  a  result 


a  =  b  = 


(7.13) 


The  dimensional  ratios  then  became 


2  2  2 

,u.  v.  ,w.  _  uv  uw  1 

(a>  *  <->  -  <->  -  —  -  ~2  -  2 

a  b 


(7.14) 


The  force-stress  and  displacement-strain  relationships  can  then  be 
simplified  and  become 


F  F  F 

=  _I_  _  _2_  _3_ 

h  200  5  °2  =  200  *  °3  °  200 


(7.15) 


'1  10  *  c2  10  ’  3  10 


The  model  stiffness  matrix  K  can  be  simplified  by  substituting  the 
relations  7.14  into  matrix  7.12.  Note  that  one  cannot  in  general  use 
conditions  2  and  3  for  isotropic  behavior  to  further  simplify  matrix  7.12. 
The  member  stiffnesses  are  related  initially  as  required  for  isotropic 
behavior;  however,  they  change  in  value  dependent  upon  the  displacements 
and  may  lot  remain  equal  to  each  other.  The  member  stiffness  functions 
are  presented  in  the  next  section. 

Substitution  of  relations  7.14  into  matrix  7.12  yields 
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2k^  +  2k^  +  bV-2 


2k, 


2k6 


[K]  = 


2k, 


2k, 


2k^  +  2k^  +  4k^ 


2k, 


2k, 


2k3  +  2k6  +  4k5 


(7.16) 


7.3  MEMBER  STIFFNESS  FUNCTIONS 


The  stiffnesses  of  the  individual  members  of  the  model  we^e  selected 
such  that  it  was  possible  to  simulate  concrete  behavior.  The  general 
shape  of  a  stress-strain  curve  for  a  compressively  loaded  axis  as  observed 
from  the  testing  program  is  shown  in  Figure  35(a). 

The  stiffness  of  a  member  is  represented  by  the  slope  to  the  stress- 
strain  curve.  The  slope  of  the  stress-strain  curve  of  Figure  35(a)  has  the 
general  shape  as  shown  in  Figure  35(b). 

The  stiffness  function  (Figure  35(b))  possesses  a  negative  region; 
hence,  it  could  not  be  used  as  a  stiffness  function  for  a  member  in  tension. 
The  negative  stiffness  region  produces  a  decrease  ir»  ;tress.  A  stress- 
strain  curve  for  concrete  in  tension  does  not  possess  a  region  of  decreasing 
stress,  but  fractures  suddenly. 

The  exterior  members  of  the  model  produce  the  lateral  displacement 
and  therefore,  for  most  loading  combinations  it  was  not  always  necessary 
nor  desirable  for  the  exterior  members  to  fail  in  compression.  The  exterior 
members  could  not  possess  a  negative  stiffness.  A  negative  stiffness  would 
result  in  a  sudden  change  in  the  direction  of  the  lateral  displacements. 

No  such  sudden  changes  in  lateral  strains  were  observed  during  the  testing 


103 


STIFFNESS  STRESS 


program;  hence,  the  stiffness  function  shown  in  Figure  35(b)  was  not  used 
for  the  members  in  tension  or  the  exterior  members  in  compression. 

Generally  stresses  are  determined  from  known  strains;  hence,  it  was 
convenient  to  take  the  member  stiffnesses  to  be  functions  of  the  strains. 
The  mathematical  expression  for  the  stiffness  functions  (Figure  35(b))  for 
the  interior  members  (1-1,  2-2,  3-3)  in  compression  are  of  the  form 


(1  +  d)k. 


dk2(l  +  (1  -  cn 


•W 


rt  1+1 /n 


(7.17) 


where 


is  the  initial  stiffness  of  the  member 
is  a  constant  which  determines  the  influence  of  the 
individual  terms  upon  the  stiffness  k 
is  a  constant  which  affects  the  stress  drop-off  portion 
of  the  stress-strain  curves  once  the  peak  stress  has 
been  reached 

are  integers  which  determine  the  rate  at  which  the 
stiffness  changes 

are  values  of  the  variables  at  which  the  stiffness  k 
changes  to  a  different  value 

g(x),  f(x)  are  functions  of  the  model  displacements.  They  are 
listed  later. 

The  general  shape  of  the  above  stiffness  function  is  shown  below. 


n,  nx 


V  X1 
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g(x) 


The  stiffness  functions  for  the  exterior  members  in  compression  and 
of  all  members  in  tension  is  of  the  form 


(7.18) 


where 

k^  is  the  initial  member  stiffness 
n  is  an  integer  influencing  the  rate  of  change  of  k 

h(x)  is  a  function  of  the  displacements 

is  the  value  of  h(x)  at  which  the  stiffness  k  changes  value. 
The  general  shape  of  this  function  is  shown  below. 
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h(x) 


The  functions  g(x)  and  h(x)  for  the  different  elements  in  compres¬ 
sion  are  given  below.  The  stiffness  functions  require  that  these  variables 
be  positive;  hence,  the  absolute  values  are  shown. 

Member  Stiffness  h(x) 


|Xj  +  x2  -  2X3 | 


k 


3 


k 


6 


and 


Member  Stiffness 


|X2  +  X3  -  2X3!  '7.19) 

|XX  +  X3  -  2X2| 

|XX  -  min(X2,  X3)| 

|X2  -  min(Xlt  X3)|  (7.20) 
|X3  -  min(X1,  X2)| 
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The  functions  h(x)  for  the  elements  in  tension  are  given  below. 


Member  Stiffness 


h(x) 

|X2|  (7.21) 

|x3| 

.  707  +  X2  +  X3| 

.707|X1  +  X2  +  X3|  (7.22) 

.707|X1  +  X2  +  X3| 


The  variables  used  in  the  above  stiffness  functions  in  compression 
were  selected  such  that  the  model  would  yield  an  indefinite  strength  for 
the  triaxial  case  of  equally  applied  loads.  This  is  in  conformity  with 
test  results. 


7.4  SOLUTION  OF  THE  MODEL  EQUATIONS 
The  model  equations  are 


F  =  ATkAX  «*  KX 


(7.23) 


The  matrix  K  is  a  function  of  the  displacements,  hence 

F  -  K(X)X  (7.24) 
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The  above  equations  are  easily  solved  provided  that  Lhe  displace¬ 
ments  X  are  given  and  the  forces  F  are  determined;  however,  in  the 
development  of  the  stiffness  values  and  control  parameters  utilized  in 
the  model,  it  was  necessary  to  determine  the  displacements  corresponding 
to  a  given  set  of  forces.  Inversion  of  the  above  equations  was  difficult 
if  not  impossible  due  to  the  nonlinear  nature  of  K(X).  This  difficulty 
was  avoided  by  using  a  procedure  suggested  by  Richard  and  Goldberg  (19). 
They  suggested  converting  nonlinear  equations  to  a  set  of  linear  ordinary 
first  order  differential  equations.  Equations  7.24  are  then  converted 
to  the  differential  equations: 

dF  =  K(X)dX  (7.25) 

Methods  exist  for  numerically  solving  these  differential  equations. 

For  simplicity,  Eulers  point-slope  method  was  used.  Equation  7.25  then 
becomes 


AF  -  K(X)AX 


(7.26) 


with 


F 


i+1 


Fj^  +  AF 


(7.27) 


X 


i+1 


X£  +  AX 
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where  the  equations  now  yield  the  load  and  displacement  increments  in¬ 
stead  of  the  total  load  and  total  displacement. 

More  accurate  and  sophisticated  methods  of  numerical  integration 
could  be  used;  however,  Eulers  method  is  the  simplest.  Using  Eulers 
method  does  restrict  the  size  of  the  increments.  They  should  be  kept 
small  to  reduce  errors  in  integration.  A  more  accurate  integration  pro¬ 
cedure  may  be  utilized. 

7.5  MODEL  STIFFNESS  CONSTANTS 

The  initial  member  stiffness  values  were  selected  such  that  the 
model  would  closely  simulate  the  behavior  of  concrete  as  observed  from 
a  uniaxial  test.  The  relationships  between  test  values  and  the  initial 
model  stiffness  values  were  obtained  utilizing  the  following  conditions: 

1.  Isotropic  behavior  of  model  (see  Section  7.2). 

2.  The  initial  modulus  of  elastiticy  Eq  and  the  initial  Poisson's 
ratio  y  are  known  from  a  uniaxial  test. 

3.  During  the  first  stress  and  strain  increment,  the  member  stiff¬ 
nesses  are  constant. 

These  conditions  are  not  sufficient  to  determine  the  entire  stress-. strain 
curves.  They  control  only  the  initial  slopes  to  the  stress-strain  curves. 
The  factors  which  determine  the  shape  of  the  entire  stress-strain  curves 
are  presented  in  Section  VIII. 

The  model  equations  7.26  can  be  written  as 
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where  are  the  elements  of  the  model  stiffness  matrix. 

For  the  uniaxial  case,  AF2  =  AF^  *  0  and  AX2  =  AX^.  The  above 
equations  then  reduce  to 
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(7.29) 


Now  Poisson's  ratio  is  defined  as 


lateral  strain 
axial  strain 


hence 


(7.30) 


Note  that  for  AX2  to  be  equal  to  AX^,  the  following  relationships 
must  be  initially  satisfied: 
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K12  =  K13 


J.31) 


K22  =  K33 


These  two  conditions  are  satisfied  by  the  isotropic  requirements 
mentioned  earlier.  Using  the  isotropic  requirements  regarding  the  ele¬ 
ment  stiffness,  the  model  matrix  7.12  becomes 


4k^  +  4k2 


2k, 


2k, 


4k^  +  4k2 


2k, 


2k, 


2k, 


2k^  4k^  +  4k2 


(7.32) 


remembering  that  this  specialized  matrix  is  valid  only  for  the  first 
stress  and  strain  increment. 

The  following  mathematical  development  also  applies  only  to  the  first 
stress  and  strain  increment. 

Recalling  Equation  7.29(b). 
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Now 


I AX2 1  ~h2 

T^T  ^22  +  ^23 


(7.33) 


Substituting  the  valuer  for  K^2,  ^22*  anc*  K23  *nto  t^ie  a^ove  equation 
yields 


+  2k, 


(7.34) 


The  second  relationship  is  obtained  by  considering  the  initial 
modulus  of  elasticity  as  obtained  from  a  uniaxial  test.  The  initial 
modulus  of  elasticity  is 


or 
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(7.35) 


Now  recall  Equation  7.29(a) 


AFX  *  KUAX1  +  (K12  +  K13)AX2 


or 
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(7.36) 


AF1  M2 
AXX  K11  +  (K12  +  K13*  AX: 


Recall  that 
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(7.37) 


Substitution  of  the  stiffnesses  for  K^,  K^»  K^»  K 22 »  and  yields 

2uE  -  4k,  +  4k0  -  2k,  y 
o  1  2  1 


or 


uEq  -  2(kx  +  k2)  -  kxy  -  (2  -  y)kx  +  2k2 


(7.38) 


Recall 


3k1  +  2k2 


Solving  Equations  7.34  and  7.38  for  k^  and  k2  yields: 
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2(1  -  u  -  2y2) 


(7.40) 


k 


The  two  above  relationships  yield  the  initial  stiffnesses  and 
k^  when  given  the  initial  modulus  of  elasticity  and  initial  Poisson 
ratio  from  a  uniaxial  test. 

Again  the  mathematical  development  in  this  section  applies  only  to 
the  first  stress  and  strain  increment. 
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SECTION  VIII 


MODEL  CHARACTERISTICS 


8.1  CONTROL  PARAMETERS 

A  number  of  control  parameters  are  incorporated  into  the  model. 

These  parameters  determine  the  peak  stress  values  and  the  strain  at  the 

i 

peak  stresses  as  well  as  the  general  shape  of  the  predicted  stress- 
strain  curves.  The  function  of  the  control  parameters  is  to  change  the 
stiffnesses  of  the  individual  members  of  the  model.  The  parameter 
values  were  determined  by  the  method  of  successive  attempts,  that  is,  for 
a  given  loading  proportion,  the  parameter  values  were  changed  until  the 
model  approximately  predicted  the  test  results.  The  model  was  required 
to  approximately  predict  the  experimental  stress-strain  curves  for  all 
loaded  axes.  Control  parameter  values  were  obtained  for  several  loading 
proportions.  An  equation  was  then  fitted  to  the  parameter  values.  This 
empirical  equation  was  then  incorporated  into  the  computer  program  for 
the  model  solution  and  again  verified  by  comparing  the  model  predicted 
results  to  the  test  results.  In  general  the  control  equations  were 
different  for  each  loading  region,  that  is,  a  set  of  control  equations 
was  used  for  the  compression-compression-compression  region  and  a  dif¬ 
ferent  set  of  control  equations  was  used  for  the  compression-compression- 
tension  region. 

The  control  equations  (listed  later  in  this  section)  are  functions  of 
the  stress  increment  ratios  Aa^/Ao^  and  AOg/Ao^.  The  stress  increments  are 
not  initially  known  for  a  given  set  of  strain  increments;  hence,  the  control 
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parameters  are  determined  from  the  previous  set  of  stress  increments.  It 
is  for  this  reason  that  several  passes  through  the  computer  program  for 
model  solution  must  be  made  to  obtain  a  better  result.  This  difficulty 
may  possibly  be  lessened  by  using  a  more  accurate  integration  procedure. 

It  would  have  been  more  desirable  to  express  the  control  equations 
as  functions  of  the  strains;  however,  this  was  not  practicable  during  the 
model  development  stage.  The  predicted  strains  changed  considerable  when¬ 
ever  the  control  constants  were  changed. 

The  control  parameters  are  listed  below  along  with  their  effect  on 
the  model  results.  The  terminology  used  in  conjunction  with  the  parameter 
descriptions  with  regard  to  a  typical  stress-strar.n  curve  is  indicated  in 
Figure  36.  The  parameter  names  are  the  same  as  used  in  the  computer 
program  for  the  model. 

X10  -  This  parameter  controls  the  stiffness  of  the  exterior  diagonal 
members  of  the  model.  It  has  a  pronounced  effect  on  the  gener¬ 
al  slope  of  the  intermediate  part  of  the  stress-strain  curve 
and  on  the  magnitude  of  the  lateral  deformations.  It  was  used 
primarily  to  regulate  the  lateral  deformations. 

X0C2  -  This  parameter  partially  determines  the  stiffness  of  the  inter¬ 
ior  members  of  the  model.  Particularly,  it  controls  the  stiff¬ 
ness  value  of  the  second  term  of  the  stiffness  function  7.17 
listed  in  Section  VII.  It  determines  the  strain  at  which  the 

-D 

peak  stress  occurs.  It  also  affects  the  magnitude  of  the  peak 
stress  and  the  range  of  the  intermediate  part  of  a  predicted 
stress-strain  curve. 
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Figure  36.  Strain  Regions  of  a  Compressive  Stress- Strain  Curve 


DCN  -  This  parameter  governs  the  relative  effect  of  the  individual 
terms  of  the  stiffness  function  7.17.  It  affects  the  slope 
of  the  intermediate  part  as  well  as  the  drop-off  of  the  stress- 
strain  curve.  The  only  stiffness  function  that  can  take  on 
negative  values  is  the  second  term  of  the  stiffness  function 
7.17. 

X21  -  This  parameter  controls  the  value  of  the  first  term  of  the 
stiffness  function  7.17.  It  determines  the  strain  at  which 
the  interior  members  decrease  in  stiffness.  The  value  of  this 
parameter  determines  the  point  at  which  the  stress-strain  curve 
changes  from  the  initial  part  to  the  intermediate  part.  It 
also  determines  the  range  of  the  initial  part  and  influences 
the  range  of  the  intermediate  part  of  the  stress-strain  curves. 

ET  -  This  parameter  governs  the  strain  at  which  the  model  members 
fail  in  tension.  It  controls  the  maximum  tensile  strength 
that  is  obtained.  Whenever  a  member  fails  in  tension,  the 
stress  in  that  member  must  be  set  to  zero.  A  zero  stiffness 
will  result  in  a  constant  stress  level  which  is  not  consistent 
with  actual  concrete  behavior. 

CON  -  This  parameter  controls  the  drop-off  part  of  the  stress-strain 
curves.  It  also  partially  controls  the  total  stress  decrease 
and  the  slope  of  the  downward  part. 

XOT  -  This  parameter  has  a  value  equal  to  or  slightly  less  than  ET. 

Its  function  is  to  provide  the  slight  nonlinear  part  of  a 
tension  stress-strain  curve  near  the  ultimate  tensile  strength. 
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The  control  parameter  equations  are  listed  below  along  with  the  loading 
regions  they  control.  The  notation 

B  =  Aa3/Aa1 

G  =  A a2/Ao1 

is  used  to  simplify  the  listing  of  the  control  equations. 

(a)  Biaxial  compression 
X10  «  0.055 

X0C2  *  0.0265  +  ('(0.156  +  (G  -  0.1)  (-0.685  +  (G  -  0.2)  (1.638 

+  (G  -  C.4)(-2.81  +  (G  -  0.6)0.705  +  (G  -  0.8)  (-4.25 

+  (G  -  0.9) (3.867)))))))) 

DCN  =  1.0  -  0.5G 
X21  «  0.03  -  0.0125G 
ET  =  0.0008 
CON  =  1.05  +  0.5G 
XOT  =  0.0008 

(b)  Triaxial  compression 

X10  =  0.05  +  0.5B  +  3. 0B  (B  -  0.1)  +  (-0.06  +  1.80B)G 
-(0.05  +  1.50B  +  7.50B  (B  -  0.1))G(G  -  0.2) 

XOC2  =  0.035  +  1.65B  +  3.30B  (B  -  0.1)  +  (0.022  -  0.72B  +  23.35B 

(B  -  0.1) )G  +  3.75B  (B  -  0.1)G(G  -  0.2)  -  50. 0B  (B  -  0.1) 
G(G  -  0. 2)  (G  -  0.4) 

DCN  =  1.0  -  6.95B  +  34.75B(B  -  0.1)  -  4.0BG  +  5.0BG(G  -  0.2) 

X21  =  (0.04  -  0.01B  -  0. 25B(B  -  0.1)  tanh(G(9.70  -  57. 0B  +  235. 0B 
(B  -  0.1))) 


i  > 
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ET  =  0.008 


CON  -1.0 
XOT  =  0.0008 

(c)  Uniaxial  compression,  tension;  Biaxial  compression-tension; 
Triaxial  compression-compression-tension,  compression-tension- 
tension 

X10  =  0.055 
X0C2  =  0.0265 
DCN  -  1.0 
X21  -  0.055 

ET  -  ((0.02671  +  0. 1668B)G)  exp((-20.383  +  7. 773B)G(1‘ 25  "  °*5B)) 
+  0.0012 
CON  =  1.0 
XOT  =  0.0008 

(d)  Biaxial  tension;  Triaxial  tension 
X10  -  0.055 

X0C2  -  0.0265 
DCN  -1.0 
X21  -  0.055 

ET  -  0.0012  -  0.0001G(3  +  3.25B) 

XOT  -  0.0008 

Two  integers,  designated  NT  and  NC  in  the  computer  program,  also  have 
an  influence  on  the  shape  of  the  predicted  stress-strain  curves.  NT  and  NC 
are  the  names  used  in  the  computer  program  to  represent  the  integers  n  and 
n^  used  in  conjunction  with  the  member  stiffness  functions  described  in 
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Section  VII.  The  values  of  these  integers  affect  the  quickness  or  sharp¬ 
ness  of  the  transition  zones  between  the  initial  and  intermediate  part  and 
in  the  vicinity  of  the  peak  stresses. 

A  detailed  description  of  the  computer  program  for  the  model  is 
presented  in  Appendix  I. 

8.2  COMPARISON  OF  MODEL  PREDICTED  RESULTS  WITH  TEST  RESULTS 

The  model  predicted  results  and  the  test  results  for  several  different 
loading  proportions  are  shown  in  Figures  37  through  50.  Stress-strain 
curves  obtained  from  tests  are  shown  along  with  the  model  predicted  re¬ 
sults  for  typical  cases.  The  comparison  between  model  predicted  and  test 
results  for  compressive  loads  are  generally  quite  good  but  the  model  pre¬ 
dicted  results  for  the  intermediate  and  minor  stresses  are  usually  slightly 
larger  than  the  test  results. 

The  comparison  between  model  and  test  results  are  generally  not 
good  for  cases  where  a  tensile  stress  is  applied  to  an  axis.  There  was 
also  a  problem  in  measuring  the  smaller  strains  associated  with  tensile 
loading.  The  extensometers  were  designed  for  large  strain  measurements 
and  may  not  be  very  accurate  for  small  strains. 

The  model  does  not  usually  produce  lateral  deformations  as  large  as 
those  observed  from  tests.  The  model  geometry  is  such  that  a  maximum 
possible  ratio  of  minor  strain  to  major  strain  is  unity.  This  does  not 
permit  the  large  lateral  deformations  that  were  observed  in  the  direction 
of  the  minor  stress. 

The  comparison  is  not  always  very  close  after  the  peak  stress  had 
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been  obtained.  The  model  continues  predicting  stresses  and  strains  due 
to  the  given  loading  proportions.  In  conducting  the  tests,  it  was  impos¬ 
sible  to  maintain  a  given  loading  proportion  once  the  peak  stresses  had 
been  obtained.  A  comparison  of  predicted  and  test  values  beyond  the  peak 
stresses  is  therefore  not  logical. 

8.3  MODEL  LIMITATIONS 

The  model  was  developed  in  principal  stress  and  strain  coordinates; 
hence,  the  model  itself  predicts  principal  stresses  from  principal  strains. 
For  coordinate  systems  other  than  principal  stress  and  strain  coordinate 
systems,  the  appropriate  stress  and  strain  transformation  relations  must 
be  used  in  conjunction  with  the  model. 

The  control  equations  were  based  on  a  limited  number  of  tests;  there¬ 
fore,  a  limited  number  of  different  loading  proportions.  In  triaxial 
compression,  failure  was  not  accomplished  for  minor  to  major  stress  ratios 
above  0.2.  The  model  predicted  results  are  therefore  not  verified  by  tests 
for  minor  to  major  stress  ratios  above  0.2.  The  accuracy  of  predicted 
stress-strain  values  obtained  by  extrapolation  beyond  the  value  of  0.2  is 
unknown. 

There  was  a  relatively  large  number  of  test  results  available  in  the 
compression  region  within  the  limitations  mentioned  above;  however,  in  the 
cases  where  there  was  tension  on  at  least  one  of  the  axes,  the  number  of 
test  results  was  limited.  Tension  tests  were  much  more  difficult  and  time 
consuming  than  were  the  compression  tests.  The  region  between  test  points 
is  greater  and  not  all  of  these  gaps  in  test  results  were  checked.  The 
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model  probably  yields  results  similar  to  results  obtained  from  tests. 

During  testing,  no  sudden  changes  in  concrete  behavior  were  detected  within 
a  given  region. 

The  testing  program  utilized  essentially  one  concrete  strength;  how¬ 
ever,  there  was  batch  to  batch  variation.  It  is  believed  that  the  batch 
to  batch  variation  is  small  in  comparison  to  the  variation  that  would  have 
been  obtained  using  different  concrete  strengths.  In  the  model  develop¬ 
ment,  the  average  result  was  used;  hence,  this  would  relate  to  an  average 
concrete  strength. 

8.4  OTHER  CONCRETE  STRENGTHS 

The  model  control  parameters  were  based  an  experimental  stress-strain 
curves  for  tests  on  a  single  nominal  concrete  strength.  For  the  model  to 
be  of  a  general  use,  it  would  be  helpful  if  it  could  be  easily  adjusted  to 
predict  stress-strain  values  for  any  concrete  strength.  The  effect  of 
different  concrete  strengths  and  the  effect  of  different  types  of  aggregate 
used  in  the  concrete  can  be  determined  by  additional  tests  only. 

From  limited  evidence  (Part  II),  it  appears  that  by  using  the  initial 
tangent  modulus  for  a  concrete  of  any  strength,  the  initial  slopes  to  the 
predicted  stress-strain  curves  will  compare  favorably  to  the  initial  slopes 
of  experimental  stress-strain  curves.  The  peak  stresses  and  the  strains 
at  the  peak  stresses  as  well  as  the  shape  of  the  predicted  stress-strain 
curves  are  determined  by  the  control  equations.  Indications  are  that  the 
control  equations  must  be  altered  if  the  model  is  to  closely  predict  stress- 
strain  values  for  different  concrete  strengths.  Adjustments  would  have  to 
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be  made  to  alter  the  peak  stress,  the  strain  at  the  peak  stress,  and  the 
shape  of  the  stress-strain  curves.  Changing  the  control  equations  is  no 
small  task.  The  nature  of  the  control  equation  alterations  for  different 
concrete  strengths  and  possibly  different  concrete  mixes  is  not  known  at 


this  time. 


O  Modal  Predicted  Results 


Figure  37.  Uniaxial  Compressive  Stress-Strain  Curves 
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Figure  38.  Biaxial  Compression  Stress-Strain  Curves 
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O  Model  Predicted  Result* 
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Figure  41.  Biaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  42.  Biaxial  Tension  Stress-Strain  Curves 
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Figure  43.  Triaxial  Compression  Stress-Strain  Curves 
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Figure  45.  Triaxial  Compression  Stress-Strain  Curves 


Figure  46.  Triaxial  Compression  Stress-Strain  Curves 
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Figure  47.  Triaxial  Compression  Stress-Strain  Curves 
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Figure  49.  Triaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  50 


Triaxial  Tension  Stress-Strain  Curves 


SECTION  IX 


CONCLUSIONS  AND  RECOMMENDATIONS 


9.1  CONCLUSIONS 

This  study  was  undertaken  to  determine  the  behavior  of  plain  concrete 
under  combined  stresses  and  to  formulate  computer  orientated  constitutive 
relations  for  concrete.  One  nominal  concrete  strength  was  tested  under 
various  loading  combinations.  The  loading  combinations  included  the  uni¬ 
axial,  biaxial,  and  triaxial  states  of  stress  with  various  combinations 
of  compressive  and  tensile  stresses.  The  test  specimens  were  three- inch 
cubes.  The  test  information  was  obtained  as  stress-strain  records  for  the 
three  principal  directions  of  the  cubical  test  specimens.  Generally  the 
loads  were  applied  such  that  the  ratio  of  the  loads  relative  to  each  other 
was  constant.  In  a  few  cases,  an  incremental  loading  was  used. 

Several  items  were  discovered  or  observed  during  the  course  of  the 
experimental  investigation.  The  effect  of  friction  in  testing  cubes 
was  considerable.  A  much  higher  apparent  strength  was  obtained  in  tests 
not  utilizing  friction  reduction  methods;  consequently,  friction  reducing 
pads  were  utilized  throughout  the  testing  program.  The  friction  reducing 
pads  consisted  of  two  polyethelene  sheets  with  grease  between  them.  These 
pads  were  placed  between  the  platens  and  cube  surfaces  in  the  compression 
tests. 

During  the  testing  program  it  was  noted  that  somewhat  different 
strength  values  were  obtained  depending  upon  the  orientation  of  the  cube 
with  respect  to  the  applied  loads  and  the  direction  of  casting.  The  cubes 
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were  not  isotropic.  Visual  examination  of  the  cubes  clearly  indicated 
that  the  orientation  of  the  coarse  aggregate  was  not  the  same  in  the 
direction  of  casting.  It  appears  that  spherical  aggregate  would  have  to 
be  used  to  eliminate  the  anisotropic  condition.  The  tensile  strength  in 
the  direction  of  casting  was  approximately  30  percent  less  than  the  ten¬ 
sile  strength  in  one  of  the  other  directions.  A  similar  observation  was 
made  by  Ash  (31)  during  a  study  in  which  specimens  were  loaded  in  direct 
tension.  The  tensile  strength  in  the  direction  of  casting  was  approxi¬ 
mately  24  percent  less  than  the  tensile  strength  in  the  direction  per¬ 
pendicular  to  casting.  He  attributed  the  difference  in  strength  to 
bleeding  in  concrete  which  produced  weak  regions  beneath  aggregate 
particles.  The  uniaxial  compressive  strength  was  not  greatly  affected 
by  the  orientation  of  the  cubes;  however,  the  biaxial  and  triaxial  test 
results  were  noticeably  affected  by  the  cube  orientation. 

The  failure  mode  of  the  cubes  was  identical  to  that  described  by 
several  investigators.  The  cubes  split,  apparently  in  tension,  in  the 
direction  of  the  lowest  compressive  stress  or  the  highest  tensile  stress. 
The  cubes  always  expanded  in  those  directions.  The  fracture  planes  were 
generally  perpendicular  to  the  direction  in  which  the  lowest  compressive 
stress  or  highest  tensile  stress  was  applied. 

It  was  noted  that  the  slightest  confinement  of  a  test  cube  in  an 
intended  uniaxial  test  would  noticeably  affect  the  maximum  strength.  Also 
in  a  biaxial  test,  the  results  would  be  affected  whenever  a  slight  lateral 
pressure  accidently  occurred  on  the  "free  surfaces."  Once  discovered. 
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measures  were  taken  to  prevent  lateral  pressure  from  developing  on  the 
unconfined  sides. 

In  triaxial  compressive  tests,  it  was  noted  that  the  maximum  stress 
was  greatly  affected  by  the  magnitude  of  the  minor  stress  and  somewhat 
affected  by  the  intermediate  stress.  The  strains  were  generally  affected 
in  the  same  manner.  The  largest  minor  stress  that  could  be  applied  and 
still  obtain  failure  was  approximately  7000  psi.  A  minor  stress  above  this 
value  would  prevent  failure  of  the  test  specimen  before  the  capacity  of  the 
testing  machine  was  reached. 

Whenever  applicable,  the  test  results  from  this  investigation  were 
compared  to  the  test  results  reported  by  other  investigators.  The  scatter 
in  results  reported  is  rather  large;  however,  many  of  the  investigators 
used  different  testing  procedures,  equipment,  and  different  shaped  test 
specimens.  The  results  of  this  investigation  were  bounded  by  the  results 
reported  by  other  investigators. 

A  model  to  predict  concrete  behavior  was  developed.  The  mathematical 
development  of  the  model  was  theoretical;  however,  empirical  results  were 
incorporated  into  the  model  such  that  the  test  results  were  simulated. 

The  model  was  used  to  predict  the  constitutive  relations  for  concrete 
subjected  to  combined  loads.  The  loads  can  be  tensile  or  compressive. 

The  model  solution  was  programmed  for  use  with  a  computer  and  written  in 
Fortran  IV  language. 

The  model  was  developed  using  test  results  for  a  concrete  strength  of 
4000  psi.  The  test  results  in  the  compression  range  were  limited  by  the 
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minor  stress  level.  The  model  is  controlled  by  parameters  which  control 
the  shape  of  the  stress-strain  curves  and  the  maximum  stresses. 

9.2  RECOMMENDATIONS 

The  model  itself  can  be  made  to  more  accurately  simulate  concrete 
behavior.  The  degree  of  simulation  becomes  a  matter  of  refinement  in 
the  control  parameters.  The  area  of  doubt  with  respect  to  test  results 
was  in  the  tension  region.  The  accuracy  of  the  extensometers  should 
be  improved  for  this  region  and  a  testing  procedure  and  machine  devised 
such  that  bending  stresses  are  not  induced  into  the  test  specimen  during 
testing. 

The  model  was  developed  using  one  concrete  strength  only;  hence,  the 
model  predicted  results  are  applicable  to  the  one  strength  only.  The 
model  should  be  verified  for  different  concrete  strengths.  This  can  be 
done  only  through  a  testing  program. 

The  effect  of  creep  and  loading  rate  on  concrete  under  combined 
stresses  is  not  yet  established.  The  effect  of  loading  and  unloading 
concrete  under  combined  stresses  has  not  been  established  either.  The 
constitutive  relations  will  not  be  completely  general  until  all  of  these 
effects  are  known. 
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APPENDIX  I 


COMPUTER  PROGRAM  FOR  THE  MODEL 


1.1  GENERAL 

The  computer  program  for  the  model  which  was  developed  and  described 
in  Sections  VII  and  VIII  is  presented  in  detail  in  this  appendix.  Fortran 
VI  was  used  in  the  progra  ming.  The  order  of  presentation  is  as  follows: 

(a)  Driver  Program 

(b)  Subroutine  MATER6 

(c)  Sample  Input 

(d)  Sample  Output 

1.2  DRIVER  PROGRAM 

The  purpose  of  the  driver  program  is  to  read  or  compute  the  informa¬ 
tion  that  is  transferred  into  subroutine  MATER6.  In  addition,  the  driver 
program  prints  or  writes  information  computed  within  the  subroutine. 

Entry  into  the  subroutine  is  made  through  a  CALL  statement.  Note  that  the 
driver  program  as  presented  is  self-contained  and  that  it  must  be  modified 
for  use  in  a  more  general  program.  A  data  card  with  a  strain  of  unity 
inserted  at  the  end  of  a  set  of  strain  data  cards  will  reset  the  program 
to  accept  another  set  of  strain  data  cards.  A  data  card  with  strains  larger 
than  unity  will  terminate  the  program.  The  driver  program  also  equates 
new  stresses  and  strains  to  old  stresses  and  strains  after  a  set  of  new 
stresses  have  been  computed. 

Included  in  this  section  is  a  list  of  the  variables  used  in  the 
driver  program,  a  flow  chart,  and  a  listing  of  the  driver  program. 
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DRIVER  PROGRAM 


BU 

E 

SK(I) 
SIG1(I) 
SIG(I) 
EPS1(I) 
EPS (I) 


LIST  OF  VARIABLES  AND  THEIR  DEFINITION 

Initial  Poissons'  Ratio 

Initial  Tangent  Modulus 

Initial  Stiffness  of  Model  Members 

Old  Stresses 

New  Stresses 

Old  Strains 

New  Strains 
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FLOWCHART  FOR  DRIVER  PROGRAM 
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002*5  SIGI  (J)  =  S1G(J) 

0026  920  FPill J)=FPSI J) 

0027  Go  TO  105 

0028  FNn 


1.3  SUBROUTINE  MATER6 


The  purpose  of  Subroutine  MATER6  is  to  compute  a  set  of  stresses  from 
a  given  set  of  strains.  The  subroutine  utilizes  incremental  principal 
strains  from  which  incremental  principal  stresses  are  determined.  Con¬ 
verting  to  principal  strains  and  stresses  is  a  necessity  and  must  be  ac¬ 
complished  before  strain  increments  can  be  computed.  Principal  strains 
are  determined  through  the  use  of  the  standard  subroutine  EIGEN.  Sub¬ 
routine  EIGEN  is  also  used  to  determinj  the  directional  cosines  of  the 
principal  new  strains.  Subroutine  EIGEN  also  sorts  the  principal  stresses 
or  strains  from  largest  to  smallest.  The  principal  old  stresses  are 
determined  using  stress  transformation  equations  and  the  directional  cosines 
obtained  from  the  old  strains.  The  principal  incremental  stresses  are 
determined  and  Subroutine  MATER6  computes  the  incremental  stresses  and 
total  stresses  and  then  converts  the  stresses  back  to  the  original  co¬ 
ordinate  system.  Control  is  then  transferred  back  to  the  driver  program. 

Caution  must  be  exercised  in  the  selection  of  the  strain  increments. 
Subroutine  MATER6  utilizes  Euler's  point-slope  integration  method;  hence, 
it  will  become  unstable  if  the  strain  increments  are  not  kept  sufficiently 
small.  The  program  does  contain  an  accuracy  check,  and  if  exceeded, 
divides  the  given  strain  increment.  This  should  not  be  considered  to 
completely  safeguard  against  an  unstable  solution. 

The  list  of  variables  used  in  Subroutine  MATER6  is  included  in  this 
section  along  with  the  flowchart  and  the  computer  program  listing. 
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SUBROUTINE  MATER6  LIST  OF  VARIABLES 


A(I) 

AK(I) 
ASAT (I, J) 
B 

CK 

CON 

DC(I) 

DCN 

DCP 

DDF (I) 

DF(I) 

DX(I) 

EIGEN 

G 

G1,G2,G3, 

G4 

IP 

NC,NT 

NCR 

NTIME.IQ 

01,02,03 

PSNO(I) 

PSNN(I) 


Temporary  v'lues  of  stresses  and  strains 
Member  stiffness  values 
Model  stiffness  elements 
Ratio  of  DDF (3) /DDF  (1) 

Temporary  stiffness  value 
A  control  parameter 
Directional  cosines 
A  control  parameter 

A  fraction  of  the  stiffness  of  the  compressive  members 

Sorted  absolute  values  of  DF(I),  largest  to  smallest 

Principal  incremental  stresses 

Model  incremental  displacements 

Subroutine  for  computing  eigenvalues 

Ratio  of  DDF (2) /DDF (1) 

Control  equation  variables 
Counters 

Integers  used  in  stiffness  functions 
Counter  determining  the  stress  condition 

Counter  determining  the  number  of  passes  through  model  computations 

Control  equation  variables 

Principal  strains  obtained  from  old  strains 

Principal  strains  obtained  from  new  strains 
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R(I) 

SIGP(I) 

SIGS(I) 

TEF 

V 

xoc 

XOC2 

XOT 

X21 

X(I) 


Temporary  values  of  directional  cosines 
Principal  stresses  from "old  streSses  ' 

Computed  principal  stresses 
Temporary  principal  stress 

Displacements  obtained  from  new  principal  strains 

A  control  parameter 

A  control  parameter 

A  control  parameter 

A  control  parameter 

Model  displacements  from  old  strains 


FLOW  CHART  FOR  SUBROUTINE  MATER6 
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FLOW  CHART  FOR  SUBROUTINE  MATER6  (cont'd) 


153 


FLOW  CHART  FOR 


Check  signs 

of  DF(I) . 

Increase  counter  NCR 

as  follows: 

OF  (I) 

NCR 

Pos. 

1 

Zero 

2 

Neg. 

3 

Initial  values  of 
control  constants 
DCN ,  CON,  XOC2 , 
X21 ,  X10,  ET 


FLOW  CHART  FOR  SUBROUTINE  MATER6,  (cont'd) 


Reevaluate  control 
constants  for: 
Uniaxial  compression 
Uniaxial  tension 
Biaxial  T-C 
Triaxial  T-C-C 
Triaxial  T-T-C 


ET 


i 

»  i 
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FLOWCHART  FOR  DETERM  I  NAT  TON 
OF  MEMBER  STIFFNESSES 


i 
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FLOW  CHART  FOR  DETERMINATION  OF 
MEMBER  STIFFNESSES,  (contd) 
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SUBROUTINE  MATER 6  LISTING 
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SUBROUTINE  MATER 6  LISTING  (cont’d) 
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SUBROUTINE  MATER6  LISTING  (cont'd) 
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01O1  FNO 


1.4  SAMPLE  INPUT 


Input  to  the  computer  program  consists  of  the  following: 

(a)  Poissons  Ratio  and  Initial  Tangent  Modulus  -  BU,E 

Both  of  these  values  are  included  on  the  same  data  card 
shown  below.  Their  FORMAT  specification  is  F10.0.  Poissons 
Rati.^  's  dimensionless  and  the  unit  for  the  tangent  modulus  is 
psi.  This  is  the  first  data  card  in  the  program. 


L'.  15  3282U0U.U 

I  I 

■  |  0  0  0  0  0  0  0  0  C  0  0  0  0 III 0  |0  C  0  D  0  0  G  0  0  0  0  0  0  0  0  0  0  0  D  I)  0  0  0  0  0  0  0  0  0  Ii  0  0  0  0  0  C  0  0  0  0  0  5  0  0  0  0  I 

It*  4  &  4.  711  IQ  I.  12  I*  MISS?  f  UMZi  It  ii  A  *  71  *17  13  MM  31  iff  IS  «•«>  UM  «  U  .»  w  >i  W  >J  X  %  1#  W  i*  *  *«  w  it  u>  »• 

1 1 1|  1 1 1  mi  1 11111111  in  i  n  i  ii  n  i  n  I  n  n  l  in  i  ii  1 1  it  in  1 1 1  n  ii  i  li  i  1 1 

1222 22222222|l|22l222222222222222222222222222222Z2222 2222222222222; 
S3|3 3 3223331313333 111113333233333133333333333333323 333333333333333: 
44444  4<4444444444444444  44444444444<U4<M4  1  U  4  4  4  U  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  1 

5  5  5  5  1  5  5  5  S  5  5  S  S  5  S  5  5  5  5  5  S  5  S  5  5  5  5  5  5  S  5  5  4  1  5  5  5  i  S  b  5 1  4  5  5  5  5  5  3  5  5  5  5  S  4  1  5 1  4  5  5  5 1 3 •> ' 
6SGE6CCB6C666SKfiSSCSIESSSE6BC66CCCE6B666C66E6SB6B666BB866E6666(.bG(;i 
|  777777777717777777777777777777777777777777777777777777777777777777; 

|  l E| B £  8 E E Et 8 l  l|liS!|llll  1  3  8  8  1  1  3  8  S I  8 8  I  8 i i 8  » 8 3  S3  i  »  3  3  5  8  B 8  8 8 3 U  3  6  8  8  8  8  « 8  I 

3  9  5  J  ?  :  S  5  9  9  9  ?  9  9  9  9  9  9  9  8  3  9  9  9  9  5  9  3  9  9  3  9  9  9  9  9  9  9  9  3  9  3  9  9  9  9  3  3  3  3  9  9  3  9  3  3  9  ?  9  9  3  9  9  3  9 

I  4  j  •>  .  4  4  \  4  ii  I  :•-»»»■  li  C  aEM  »  aa»«  IS  J>  J,  «  JXK.  C  U  V  V,  i-  . 
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(b)  Old  Stresses  and  Old  Strains  -  SIG1,  EPS1 


These  values  are  read  only  at  the  beginning  of  a 
series  of  given  strains  and  are  all  zero.  Their  FORMAT 
specification  is  F10.9.  The  unit  of  the  strains  is  in/in 
and  the  unit  of  the  stresses  is  psi.  A  sample  card  is 
shown  below. 


I  I  I  I  I  I 

|o|oosoooo|0|ooooooo|o|ooooooQ|o|ooooooo|o|OQooooo|o|ooooooooooo30fi 

i  i  i  t  !  I  >  I  i  it  n  u  n  m  «  K  i)  II  rt  v  unit  an  n  n  it  n  »  js*  J>  nnuniMiMnuti  *n»wii  »  « at  »  M  97  si  j»  to  tin  dm  mm 

1  1  1  1  1  1  I  1  1  1  I  1  1 1  I  1  1  1  1  I  1  1  1  1  1  1  1  1  1  t  1  1  !  1  1  I  1  1  I  I  I  1  1  1  I  1  I  1  1  1  1  1  I  I  1  1  1  1  I  1  1  1  1  I  1  I  1 

2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  2  7  2  2  2  2  2  2?  2222222222? 

3|333333333|333333333|333333333|333333333|3^3333333|333333333333333 

4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  «  4  4  4  «  4  4  «  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4  4 

5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  S  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  i  5  5  5  S  5  5  5  S  5  5  5  5  5  5  5  5  5  5  i  5  5  5  S  1  5  5  '» !.  S  5  i  i  i 
£  t  E  6  G  5  (  6  8  6  E  6  (  S  £  6  S  S  S  6  B  t  S  6  6  6  I  6  E  6  6  S  (  S  6  E  S  6  i  I  6  E  S  E  6  6  6  6  6  6  6  E  t  fi  6  6  6  f.  b  6  6  6  fi  S  fi  fi  fi 

|  mmmnmnm  mi  mini  mini  mnmimnnn  11111  n  11111 

9999999999999999999999999999999999999999999999999999999999999999999 

1  7  3  4  5  I  1  •  •  10  11  17  13  M  IS  tt  12  tl  11  20  71  31  23  24  25  ?§  77  *1  2t  JO  39  3?  33  34  35  *  3/  31  33  40  *1  42  43  44  45  45  47  4«  41  50  51  52  53  M  55  55  57  51  5!  (0  61  U  CJ  ' 4  ’5  61  67 
IM» 
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► 


(c)  New  Strains  -  F.PS 


New  strains  are  read  in  a  sequence.  The  number  of  cards 
is  variable  depending  upon  the  number  of  strain  readings  that 
are  available  or  selected.  The  FORMAT  specification  for  strains 
and  stresses  is  F11.6  and  F11.0,  respectively.  The  units  are 
the  same  as  for  the  strains  in  (b)  above. 

The  general  strain  and  stress  tensors  are,  respectively, 
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The  corresponding  order  of  the  input  for  both  stresses  and 
strains  is 
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15  4 
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4  6  3 


163 


A  sample  input  card  is  shewn  below.  Note  that  the 

quantities  Y  ,  Y  ,  and  Y  (not  the  quantities  4  Y  » 
xy  xz  yz  z  xy 

J  Yxz,  and  ~  YyZ)  are  used  on  the  input  cards. 


.  t ! t i i j 7 1 1 1  _.  111511c'  .  luj'r’iiY.f-  -.  Ui.ihc.  iJl  .  >.•  ,i.i 

I  I  I  I  I  I 

I  I 

IO|||00|GGCO|GOOOCOOOO||0|OOROOC||00||  0O00|00iJ0fl8  0  0  0  |0OOO00D 

i  »  s  «  5  i  j  >  i  it  i.  u  n  n  a  h  a  ■  ■  »aan!»n  i«  n  Jkitx  #v  v  »*  nV>  io  r  i-  uu  ii  minx  ism  v  v  ■ 

ii  ii  i  i|i  i  in  ill  II  mm  i  n  i  m  1 1 1  ii  i  m  i  i  ii  1 1 1 1 1 1  m  (  i  1 1 1 1  n 
22227222222222222I22222222222222222I2222222222222.22222222222 

l|333333333|33333333l|3333|3333|333333333|333333333|33333333 
4444444444444444444444444444444444444M4444444444444444444n 
5  5  5  5  5  5  5  5  5  5  5  5  5  515  5  5  5  5  133  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5  5!!', 
S6CG666EG6E6868SS6StlCSSB66K^6S666|666B66R6SS66S666fi666666Ct 
7777 7l777177777777777777|777 7 7777777777777777777777777777777 
l|lllllllll|llltlllll|m  8  8  8  8  8  S|l  8  8  8  8  8  1  3  6  |  3  3  8  IS8l9l|llSM8n 

I  9  9  3  9  8  3  9  S  5  9  9  9  9  9  9  S  9  IS  !  9  9  9  9  S  S  9  9  9  3  9  9  9  3  9  9  9  9  9  9  3  S  3  3  3  3  V)  3  9  3  9  9  9  3  3  5  3 

1  a  1  4  •  •  1  g  >  wb  1:  n  m 


|  Reproduced  from 

koct  Available  copy.. 
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1.5  SAMPLE  OUTPUT 


A  sample  of  the  output  information,  which  is  the  new  stresses  SIG, 
is  shown  below.  The  strains,  which  are  indicated  as  STRAINS  (INPUT), 
produce  the  stresses  which  are  indicated  as  STRESSES  (OUTPUT).  The  strain 
values  shown  below  may  be  used  as  check  values  to  determine  if  the  computer 
program  is  functioning  properly.  The  given  strains  must  produce  the  cor¬ 
responding  stresses.  The  unit  of  strains  is  in/in  and  the  unit  of  stresses 
is  psi.  EPS(5),  SIG(5),  E?S(6),  and  SIG(6)  are  zero  and  are  not  always 


shown. 


PRINCIPAL  STRESS  RATIOS 


01  i  ^2  > 

*  l;l;l 

F  PS  ( 1 1 
SIG(l) 

FPS ( 2 ) 

sir,  (2) 

CPS(3> 
s  i  r,  1 3 ) 

FPS ( 4 ) 

s  I  r,  ( 4 ) 

S  TP  A  f  NS 
STWFSSFS 

(INPUT) 

( OUTPUT ) 

0.000010 

47. 

0.000010 

47. 

0.0000)0 

47. 

0.000000 

0. 

STRAINS 

STRFSSFS 

(  INPUT) 
(OUTPUT) 

0.000021 

98. 

0.000021 

98. 

0.000021 

98. 

0.000000 

0. 

STRAINS 
STR  eSSFS 

( INPUT) 
(OUTPUT) 

0.000034 

159. 

0.000034 

159. 

0.000034 

159. 

0.000000 

0. 

STRAINS 

STPFSSFS 

(  IN°IJT  ) 
(OUTPUT) 

0.000054 

240. 

0.000054 

240. 

0.000054 

240. 

0.000000 

0. 

STRAINS 

STPFSSFS 

( INPUT ) 
(OUTPUT) 

0.000065 

0. 

0.000065 

0. 

0.000065 

0. 

0.000000 

0. 

STRAINS 
STR  F  SSFS 

(INPUT ) 
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PRINCIPAL  STRESS  RATIOS 
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APPENDIX  II 

EXPERIMENTAL  DATA 

Included  in  this  appendix  are  the  experimental  stress-strain  curves 
that  were  not  specifically  presented  elsewhere  in  this  report.  The 
stress-strain  curves  are  shown  in  Figures  51  through  67.  Also  included 
in  this  appendix  are  Tables  8  through  10.  These  tables  contain  information 
used  in  Section  V. 


177 


c-,-0.60-. 


Figure  51.  Biaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  52.  Biaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  53.  Biaxial  Compression- Tension  Stress-Strain  Curves 
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Figure  54.  Biaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  55.  Biaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  56.  Biaxial  Compression-Tension  Strens-Strain  Curves 
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Figure  59.  Triaxial  Compression  Stress-Strain  Curves 
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Figure  61.  Triaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  62.  Triaxial  Compression-Tension  Stress-Strain  Curves 


189 


0  0.0005  Compression 

STRAIN  (in/in) 


STRAIN  {in /in) 


STRAIN  (in/in) 

Figure  63.  Triaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  64.  Triaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  65.  Triaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  66.  Triaxial  Compression-Tension  Stress-Strain  Curves 
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Figure  67.  Triaxial  Compression-Tension  Stress-Strain  Curves 
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TABLE  8 


BIAXIAL  COMPRESSION  STRENGTH  DATA 


°l/or 

°2/ar 

Nominal  Ratio 

Average 

Va2 

o,/o 

1  r 

aja 

1.20 

1.17 

1.14 

1.13 

1.03 

1.03 

1.07 

1.05 
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1.05 

10/10 

1.15 

1.14 
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1.23 

1.26 
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10/8 
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1.23 
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1.11 
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0.93 

1.24 

0.94 

1.29 

0.81 

\ 

1.18 

0.71 

10/6 

1.29 

0.78 

1.20 

0.72 
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TABLE  8  (Continued) 


°l/0r 

V°r 

Nominal  Ratio 

Average 

°l/02 

0,/o 

1  r 

V°r 

1.51 
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1.31 
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1.29 
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TABLE  8  (Continued) 


al/0r 

°2lar 

Nominal  Ratio 
°l/02 

Average 

°l/or 

o2/o. 

1.33 

0.12 

1.18 
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1.32 
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1.31 
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1.40 

0.16 

1.18 
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1.27 
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1.14 
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1.20 

0.07 
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TRIAXIAL  COMPRESSION  DATA 
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TABLE  9  (Continued) 


o 


n 


b 


to 

a 


cm 

b 


CO 

a 


to 

a 


o 

■H 

U 

jal 

<s 

ci 

b 

H 

§ 

i 

CM 

b 

•H 

i 

a 

H 

o 

$5 

b 

On 

rH 

CO 

o 

o 

o 

O 

o 

vO 

<r 

o 

vO 

o 

ON 

o 

CM 

o 

m 

r* 

r- 

H 

iH 

O 

o 

rH 

o 

O 

o 

ON 

ON 

o 

ON 

o 

ON 

o 

O 

rH 

rH 

ON 

On 

o 

ON 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

rH 

rH 

CM 

rH 

CM 

rH 

CM 

CM 

CM 

CM 

rH 

rH 

CM 

rH 

•  •  •  • 

O  CM  O'  O' 
CO  CO  O  O' 
CM  CM  CM  H 

•  •  • 

m  rH  rH 
H  H  H 
CO  CO  CO 

•  •  •  • 

O  O  O  CM 
O  O'  ©  O 
i/i  ui  m 

•  •  •  • 

<r  oo  n  co 

0  .0  O  rH 

oo  oo  cc  on 

n  cm  n  w  n  o  o 

r-»  r.  r-.  r.  co  O  O 

O'  O'  O'  O'  O'  o  o 

O  O  O  O 

o  o  o 

dodo 

o  o  o  o 

O  O  O  O  ©  rH  rH 

VO  ON  rH  vO 
00  in  O  rH 

rH  rH  CM  CO 

o  o  o 
in  m  m 

CO  CO  CO 

rH  00  O  C* 
vO  H  f^*  CM 
^  N  CO  N 

H  rH  H  M 
vo  CO  co  ON 

^  n 

CM  CM  O  cn  vD  vo  cl 

o  tn  <t  oo  oo  m 

»H  rH  CM  rH  O  CM  H 

rH  rH  rH  rH 

rH  rH  rH 

rH  rH  rH  rH 

rH  rH  rH  rH 

rH  rH  rH  rH  rH  rH  rH 

fO  O  ro 

N  CO  H 

CM  CM  CM  CO 

i/lOO 

CM  O  O 

O  H  H 

CO  sfr  rH 

in  *4*  co  co 
vo  o  in  co 

tr  vj-  m 
O'  00  00  00 
o\  r.  r.  <r 

CO  O  vo  <*  O'  o  o 
O  <r  vo  rH  CM  o  o 
n  n  in  sm  cm  oo 

rH  rH  rH  rH 

CM  CM  CM 

CO  CO  CO 

n  n  n  n 

m  tom  m  m  vo  m 

CO  O'  OO  o 
rt  00 

s*  00  t/1 

rH  rH  H 

cm  m  m 
r-* 

v£>  O  CM  I'' 

o  h  io  tn 
m  cm  O  iO 

00  ^  ^  CM 

uo  n  oo  c*. 

M1  H  H  On 

cn  cm  r-  -o  O  O 

O'  o  m  rH  o  O 

^UMOintMOO 

m  in  m  to 

vO  vO  vO 

r-.  vo  fN-  oo 

n  n  n  n  n  id  n 

o  o  o  o 

'OsHin 

OtOOOH 

O  O  O 

O  O  O 

U1  U1  U1 

O  O  O  o 

O  O  O  in 
won't 

o  o  o  o 

O  O  O  sj- 
oo  r--  n  o 

o  o  o  o  o  o  o 

00  00  rH  O  O  O  O 

vo  cn  hj-  o  oo  m  o 

tltl-t  «f 

M  M 

'j'f  >»  n 

m  m 't  in 

n  n  m-  M-  tn  M  M 

o  o  o  o 
m  m  cm  m 

rH  CM  rH  H 

O  O  O 
moo 
r-.  O  O 

O  O  O  o 

O  O  O  o 

O  ©  vo  r> 

o  o  o  o 
o  o  o  o 
r.  o  o  m 

o  o  o  o  o  o  o 
cn  o  o\  o  o  o  o 
r>»  ui  in  O  O  r*  ci 

vo  rv 

CM  O  rH  CO 
rH  rH  rH  rH 

ON  ON  On  O 
rH  rH  rH  CM 

N  N  00  O  ON  H  O 
t-i  04  04  04 

o  o  o  o 

00  o  o  o 
o  n  s  oo 

O  O  o 

O  O  O 
•»  «i  m 

o  o  o  o 

O  O  O  o 
o  <r  cm  oi 

o  o  o  o 
o  o  o  o 
m  m  vo  cm 

o  o  o  o  o  o  o 

00  VO  v£>  O  O  O  O 
Hoonnsn 

00  00  ON  o 

rH  rH  rH  CM 

rH  CM  CM 

CM  CM  CM 

O  co  c* 

CM  CM  CM  CM 

^  co  co  in 

CM  CM  CM  CM 

00  00  O'  O  O'  rH  O 
rH  rH  r-l  CM  rH  CM  CM 

/"s 

CM 

CM 

1 

CM 

i 

CM 

CM 

i 

CM 

• 

cn 

in 

i 

00 

rH 

1 

o 

i 

o 

1 

o 

d 

1 

O 

rH 

rH 

rH 

rH 

rH 

v-x 

w 

v^ 

'w' 

199 


Nominal  Ratio 


D  Oi  O  in  rn  N  H  O  n  -}  51  O  St  I''  O  O  5\5i 

\  os  ©  h  os  o  o  o'  o  h  oo  m  m  m  m  st  st 

CO  H  N  N  H  (M  N  OHrlO  OOO  Hi-IHH 

D  .  . 


O  rl  rl  N  O  O  in  rH  O  CO  O  O  Os  O  ©  O 

OiOSMSOO  H  M  O  -J  O  O  00  000(0 

Os  Os  O'  OS  O  O  •OfO'O'O  O  O  Os  OOOOs 


OOOOHH  OOOO  H  rl  O  rHHHO 


m  w  i^i  os  oo  os  cm  oo  m  m/i  in  Hriooi 

ioo«  oi-in  co  m  o  -o  o  h  h  m  n-  <n  st 

rl  Sf  H  O  O  O  N  M  n  (S|  H  rl  H  10  SO  m  U1 


HrIHOHH  OOOO  OOO  OOOO 


n 

o  o 

o 

o 

o 

o 

st 

00 

o 

rH 

m 

vO 

00  O  rH 

cn 

D 

00  SO 

in 

CO 

CN 

CO 

CN 

in 

CO 

o 

o  oo  m 

<r 

CN 

D 

r->  ov 

cn 

rH 

cn 

cn 

CN 

o 

CN 

rH 

rH 

o 

cn  rv  m 

vo 

m  i D 

m 

m 

in 

in 

rH 

rH 

rH 

rH 

CN 

CN 

CN 

st  st  cn 

m 

T> 

<u 

M 

O  CO  o  o  o 

O 

O 

00 

o 

rH 

cn  ov 

00  O  1-1  OS 

o 

D 

-J  (SI  rl  CO  O' 

00 

cn 

ON 

m 

m 

cn  cn 

o  os  m  oo 

c 

00  O  m  CM  rH 

cn 

rH 

CN 

sO 

r^. 

rH 

H  O 

CN  Mn  VO 

•H 

rH 

•  •  •  •  • 

• 

• 

• 

• 

• 

• 

•  • 

•  •  •  • 

±J 

c 

o 

m  Mn  m  m 

m 

cn 

cn 

CN 

CN 

CN 

CN  CN 

-a-  ^  cn  cn 

o 

a 

v 

OS 

/~S 

m  O  O  m  O 

o 

o 

o 

o 

o 

o 

©  o 

OOOO 

•H 

Mnmioo 

00 

so 

rH 

cn 

OS 

st  o  o 

O  SO  O  00 

W 

CO  to 

o  <t  cn  io  m 

sO 

OS 

rH 

o 

00 

CO 

in  in 

o  ^  cn  cn 

rJ 

m 

2 

d  a 

w 

st  st  st  St  st  st 

a 

rH 

A 

rH 

rs  *  *  * 

CN  CN  CN  CN 

/•N 

o  o 

o 

OOO 

o  o 

o 

CN 

o 

o 

o 

o 

OOO 

•H 

m  o 

o 

OOO 

rH  CN 

o 

m 

o 

in 

o 

o 

OOO 

2 

s 

CM  o 

sO 

m  m  cn 

rH 

vD 

r3- 

00 

CN 

rv 

cn 

m  ^  oo 

0 

p 

a  a 

a 

•IMA 

A  A 

A 

A 

A 

A 

A 

A 

AAA 

s-x 

O  c-l 

Os 

cn  cn  cn 

St  St 

cn 

st 

vO 

ON 

00 

cn 

Mnm 

2 

2 

rH 

CN  CN  CN 

rH 

rH  rH  rH 

OOOOOO  OOOO  OOO  OOOO 

mooooo  m  o  h  s  o  m  o  o  m  o  o 

cm  cm  os  m  cn  so  oo  o  o  oo  cn  oo  n  n-j  o 

«««•*««*  A  M  r  M  «  «  « 

O  N  O  (sl  (S  cn  O  O  OS  o  sO  Os  00  c0  C»  m  sO 

CN  CM  (si  CM  (M  CM  i—l  H  rH  rHHHH 


200 


TABLE  10 


TRIAXIAL 

TENSION-COMPRESSION  DATA 

o 

a_ 

a„ 

|o  | 
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1 0  1 

,0  | 
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1  *1 

1  1| 

-0.664 

-0.614 

+0.008 

0.012 

-0.924 

0 

+0.013 

0.014 

-0.759 

0 

+0.013 

0.017 

-0.498 

-0.498 

+0.012 

0.025 

-0.498 

+0.020 

+0.016 

0.033 

-0.465 

-0.465 

+0.026 

0.057 

-0.579 

-0.299 

+0.029 

0.050 

-0.616 

0 

+0.031 

0.050 

-0.661 
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+0.032 

0.049 

-0.465 

-0.232 

+0.042 

0.091 

-0.332 

-0.305 

+0.043 

0.130 

-0.462 

0 

+0.063 

0.136 

-0.265 

-0.265 

+0.050 

0.187 

-0.282 

-0.282 

+0.056 

0.200 

-0.133 

-0.133 

+0.027 

0.200 

-0.232 

-0.116 

+0.048 
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-0.279 
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+0.056 

0.200 

-0.289 

0 

+0.057 

0.200 
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+0.048 

0.200 
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0.396 

(-)  Compressive 
(+)  Tensile 
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